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An  ovci-vicv;  of  tho  5113',  ProjiVJn   io  offered  by  the  activity  chart  and  the 
accornpanyin^  (jraph  of  expenditures   in  Fijiirc  1.     Major  milestonec  uere: 
(a)    preparation  cuid  installation  of  inEtruinnntation,    (b)    initiation  of 
mcacurementc,    (c)    initiation   of  corfinutcri'/.ed   data  reduction,    and  (d) 
preparation  of  the  "Handbook  of  Solar  JLner^  Lxperiments."     The  {^-aph  of 
expenditure::  chowc  a  miifom   floy;  of  funds  after  the   initial   cost  of  the 
instruments.     The  project   remained  v.'itliin  the  budget  outlined  in  the 
orifjinal  propocal. 

A  map,  Fi^Tjre  2,    sho'./s  the  location  of  the  solar  cner;^  recording  stations. 
The  main  solar  enercy  data  bpse  cones  from  ti»e  20  recording  stations. 
About   5000  solar  data-dr^/s  liave  been  accuriulatcd  at  this  viritin^.     This 
data  is  represented  ;md  recorded  in  our  coniputer  file  and  consists  of 
about  300,000  numbers  representing;  Kolar  intensities.     Another  computer 
file,    containing  reduced  end   inte^^ratod  data,    contains  about  60,000 
niombers.     A   n\inrnation  of  this  data  in  terms  of  v.-cckly  and  monthly  averasoc 
is  given  in  this  report.      This  data  can  be  used  by  tho  solar  designer  to 
calculate  the  energy  output  of  jiis  particular  solar  convorison  device. 
Particular  excr.ples  of  applications  v.'ould  include  solar  v;ater  heaters, 
residential  and  commercial    solar  hco.ting  or  cooling  systems,    solar  prop 
dryers,    solar  lumber  drying  kilns,    and  industrial  process  heat. 

The.SII;Ii  Progi-ain  makes  use  of  volunteers  at  the  public  high  schools  to 

talwC  care  of  the  instrorncnts  and  to  mail  the  data  to  us   for  processing. 

There  are  about   50  volunteers  across  the  state  actively   involved  with 

this  program.     The  volmitecrs  are  listed  in  Figures   3  and   i\   along  v;ith 
their  schools  and  locations. 

The  recording  stiitions  often   locate  the  solar  energy  chart   recorder  in  a 
classroom  or  laboratory   room.      Teachers  have  reported  to  us  tha.t 
obse^■^^ations  of  the  solar  energy  records  have  been  incorporated  into  their 
classes.     Gome  classes  have  correlated  '.;eat}ier  observations  rnd  air 
pollution  v.'ith  the   solar   energy   recoi^ds. 


Ue   estimate  that   several   hundred   students  \;ill  have  had  some   contact  v;ith 
solar  energy  measurement   due  to  tlicir  participation  \/ith  the  rji:?.:  Prograii 
this  past  year,     "c  regard  the  solar  insolation  data  itself  as  the 
primary   goal   of  the  SIIi:  ?rogr:m   and  believe  that   this  data  \:x\\   accelerate 
solar  utilization   in  ?;ontrnr-..     The  involvement  of  the  stuf-cnts  and  their 
exposure  to   solar  ener:;y  phenomena  mry  also  contribute  significently,    in 
the  long  run,    to   solar  utilisation. 


The  fact.  tha.t  '.'c  arc  dci)Cndcnt  on  volunteers  rnci'xc  th.it  v;o  do  not  have  the 
authority  to  enforce  ctrict  operational  control'  over  the  data  utatioun. 
V.e  have  lost  zouc  data  duo  to  pov;er  outarcc,    i;chool   rcmodclin:;  projects, 
and  operational   error.     There  v.-crc  a  fa;  inctrjicc:j  of  da;na3cd  instiuments 
and  poor  electrical   connections  \;hich  resulted  in  a  tcnporr.ry  data  loca. 

Total  data  lost   from  all  ctationo   for  the  entire  year  avcraccc  aboiit  4  days 
per  month  per  station.     The  added  statistical  uncertainty  in  the  monthly 
averages  due  to  a  sanplin/;  of  8?,^  (rather  tliaa  100f.)    is  only  a  few  percent. 
Thus,   the  avera{;e  uncertainly  due  to  lo:-t  data  is  \;cll  vjithin  the  error 
boiid  of  the  inctininents  theniGclvcc  and  the  data  reduction.     A  more 
detailed  discussion  of  error  is  jivcn  in  Appendices  I  and  II. 

Many  of  the  stations-  have  a  nearly  perfect   record  for  ir.onitoring  the 
instruments  and  scndinc  in  data.     Tlic  i,-cneral  level  of  cooperation^  of  the 
voluntecrc     has  "been,    in  iny,  opinion,   very  hi^h.     I   thiri:  rr.ciny  of  the 
volunteers  are  aware  of  the  energy  crisis  and  recocni:>c  the  import;>ncc  of 
solar  eners/. 

Certificates  of  Arard  vjcro  prepared  to'..'ard  tlic  end  of  the  first  year's 
prof:rcjn.     These  are  beinc  sent  to  the  volunteer  teachers  cuid  their  hi^jh 
schools.     The  rccordinj  stations  received  vji  av;ard  as  shovn  in  Tirjure  5. 
These  v.'crc  signed   oy  John  Orth  and  Thomas  Jud^-o  as  •..'cll  ac  myself.     Cirnilar 
av/ards,   sic^ncd  only  "by  me,  \Jerc  sent  to  the  30  manual  stations.     Initial 
reports  indicate  the  success  of  this  activity. 

The  solar  insolation  data  in  the  next  section  is  the  main  body  of  this 
repoi't  and  represents  the  primarji'-  output  of  the  ?AiZ:  Procra-n.     Tccl-inical 
discussion  of  the  instruments  and  the  data  reduction  are  given  in 
Appendices  I   and  II.     Appendix  III   covers  ncv  activities  during  the 
fourth  quarter.     Appendix  IV   contains  a  lif^t  of  'U:2l  consultants. 
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SOLAR  i:nei^oy  data  n<oi:  sum' pkooijam 


Introduction 

ALL  TIIL  DATA  IN  Til?.  FOLLOVJNG  TABLP-S  TtELATTD    TO  SOUR  I1I50IJVTI0N  PALLING 
ON  A  SUIiFACE  FACHG  DIO  SOUTH  i\ND  TILTIiD  60     U?  FJ.'OM  THi:  H0RI70N.      This 
orientation  is  nearly  optimum  for  collecting  heat  energy  during  the  winter 
heating  season  in  Montajia.   . 

Total  Insolation  . 

The  tables  contain  values  of  "total"  insolation.  These  values  result 
from  the  integration  of  the  daily  solar  insolation  curves.  These  values 
represent  the  "total  energy"  that  fell  on  the  surface  of  the  transducer 
during  the  entire  day.  This  value  represents  the  mcaimum  energy  vjhich 
could  be  recovered  with  a  100f>  efficient  solar  collector.   (Real 
collectors  have  conversion  efficiencies  around  20-lCff,.) 

Threshold  Radiation 

In  addition  to  total  radiation  data,  v/e  arc  calculating  the  radiation 
above  a  series  of  thresholds.  The  reason  for  calculating  threshold  data 
derives  from  the  operating  characteristics  of  nearly  all  solar  collectors, 

A  solar  collector  is  normally  operating  only  during  periods  of  relatively 
high  solar  intensity,  due  to  the  losses  inherent  in  any  collector.  In 
early  morning  or  evening  hours  or  during  other  periods  of  lovj  radiation, 
if  the  circulation  pump  were  on,  a  collector  v;ould  lose  more  heat  than  it 
•would  gain.  The  radiation  falling  on  the  collector  must  be  of  an  intensity 
great  enough  to  raise  the  collector  surface  temperature  above  the 
temperature  of  the  storage  reservoir.  A  real  collector  will  thus  only 
capture  a  fraction  of  the  total  available  insolation. 

For  example,  if  the  storage  reseivoir  has  a  temperature  of  12^  F  and  the 
ambient  air  is  25°,  and  if  the  collector  is  dpuble-glazed,  the  heat  loss  ^ 
through  the  glazing  will  be  about  O.60  BTlj/ft^-hr-F  or  about  68  BTU/hr-ft". 
Assume  that  additional  conduction  losses  result  in  a  total  loss  of 
90  BTU/hr-ft'^.  Thus,  if  the  radiant  intensity  were  less  than  this  amount, 
the  collector  would  operate  at  a  net  LOSS. 

In  recognition  of  this  characteristic,  thc^percentage  of  total  radiation  p 
above  a  series  of  threBliolds  from  0.1  ]^' /m     to  1.2  'til /m   ,  spaced  0.1  k''/m  , 
is  tabulated.   Given  the  operating  threshold  of  a  particular  collector,  it 
is  possible  to  interpolate  the  data  to  determine  the  energy  available. 
This  value  is  further  modified  by  the  collector  efficiency  to  determine  the 
actual  energy  which  can  be  collected. 


Accuracy 

The  average  accuracy  of  the  solar  insolation  data  in  these  tables  is  -%, 
On  individual  days,    especially  overcast  days,   the  uncertainty  is  greater, 
v/hile  on  clear  days  the  uncertainty   is  letja.     Thece  uncertainties  tend  to 
average  as  the  data  itr.elf  is  averaged,     A  more  thorough  discuscion  of  the 
accuracy  of  the  instruments  and  the  data  reduction  process  is  given  in 
Appendices  I   and  II. 

llissin~  Data 

There  are  come  time  periods  in  the  taLlcs  for  particular  stationn  for 
v/hich  no  data  Mz.n  available.     These  voids  vere  due  to  a  variety  of 
causes:      instrunent  malfunction,   electric  pov;er  interruption,   chart 
paper  improperly  installed,   chart  paper  ran  out  (especially  during  the 
cummer)  I   remodeling  of  school  building,   povrcr  turned  off  during  vacations, 
data  lost  by  Post  Office,    etc.     The  tables  contain  the  statement  "14o 
Data  for  This  Period"  when  data  v.'ac  not  available. 

The  solar  desi.^ner  is  advised  to  estimate  missing  data  by  interpolating 
insolation  values  fi'om  nearby  measuring  stations.     The  map  in  Vigure  2 
v;ill  be  helpful   in  identifying  the  closest  recording  stations. 
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Solar  Insolation  Data 

Compiled  By 

Fowlkes  Enaineerin^ 

Bozemany  Montana   59715 


Locations   Billin^sr  Montana 


Interval 


DailM 
KWH/SM 


Feb  1  to  Feb  5 
Feb  6  to  Feb  12 
Feb  13  to  Feb  19 
Feb  20  to  Feb  26 
Feb  27  to  Feb  28 

Month  of  Feb 


5.41 
4.55 
3.03 
3.76 


4.20 


Mar   1  to  Mar   5   4.55 


Mar   6  to  Mar 

12 

3.94 

Mar  13  to  Mar 

19 

3.5g 

Mar  20  to  Mar 

26 

4.4| 

Mar  27  to  Mar 

31 

3.29 

Month  of  Mar 

A 

3.95 

Apt   1  to  Af»i 

Apr   3  to  API 

Apr  10  to  Apr^ 

Apr  17  to  Apr 

Apr  24  to  Apr  30   5.26 

Month  of  Apr 


X   To 


Above  Threshold  <Ky/SM) 
>,6  >.7  >.8  >.9  >1 


79  71 
64  56 
50  45 
49  43 

80  70 


06 

40 
30 

33 

56 


33 
23 
15 
24 
33 


0 
1 
0 
2 
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0- 
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82 

75 
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55 

37 

31 

10 

2 

0 

60 

60 

52 

43 

34 

13 

7 

1 

0 

69 

60 

51 

43 

33 

20 

7 

0 

0 

72 

66 

58 

49 

36 

23 

7 

2 

0 

36 

34 

32 

28 

24 

16 

0 

0 

0 

73 

53 

30 

22 

14 

8 

0 

0 

0 

0 

0 

88 

83 

77 

71 

65 

57 

46 

26 

4 

1 

0 

90 

82 

76 

69 

63 

54 

41 

13  ' 

4 

1 

0 

77 

67 

61 

55 

49 

40 

27 

9 

3 

1 

0 

94 

90 

83 

75 

66 

53 

28 

0 

0 

0 

0 

86   79   71  61 


;&      48   33   11 


May  1  to  May  7 
May  8  to  May  14 
May  15  to  May  21 
May  22  to  May  28 
May  29  to  May  31 


4.13 

95 

86 

74 

66 

58 

50 

36 

24 

10 

0 

0 

0 

4.92 

95 

92 

84 

77 

68 

60 

47 

20 

0 

0 

0 

0 

2.80 

80 

61 

50 

45 

39 

32 

25 

16 

7 

1 

0 

0 

3.96 

94 

86 

76 

67 

58 

51 

39 

15 

2 

0 

0 

0 

4.97 

96 

93 

87 

84 

72 

61 

39 

9 

0 

0 

0 

0 

Month  of  May 


4.o; 


9: 


66   57   49   37   18 


0 
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SOLAR  mraCY  DATA  FROK  IIAMUAL  STATIOIIS 


Introduction 

Thirty  hifjh  schoolc  v;cre  instrumented  v<'ith  SIMII-IM  manual  solar  energy 
instrmnents.     The  schools  were  asked  to  take  at   least  one  solar  energy- 
reading  each  school  day,    record  it  on  a  prepared  chart,    and  send  it  to . 
us  each  month  for  compilation.     The  instnimcnts  were  distributed  last 
Hay  with  the  undcrstajidintj  that  data  would  be  collected  tstarting  with 
the  beginnin.'^  of  school   in  Geptember^     Hy;:h.-iKi  and  Scobey  decided  not  to 
participate  and  the  instruments  will  be  reassigned. 

The  manual  stations  were  instrumented  with  the  assumption  that  the  data 
would  be  considerably  less  reliable  than  data  from  the  ?0  recording 
stations.     This  is  because  the  manual  data  is  taken  periodically  rather 
than  continuously  and  due  to  reading  and  orientation   errors.     Often 
different  teachers  and  students  take  the  data  during  the  month. 

Pata 

An  example  of  the  data  sheets  is  shown  in  Figure  6.  This  is  one  of  the 
more  complete  data  sets  we  have  received.  To  date  about  25-30f.  of  the 
stations  have  taken  several  readings  during  each  day,  about  25-35^.  have 
taken  more  than  one  reading,  and  about  ^?^.   have  sent  us  only  noon 
readings.  V'e  have  not  yet  received  data  from  1&,j  of  the  stations.  The 
quality  of  the  data  seems  to  be  improving  as  the  personnel  become  more 
familiar  v/ith  the  instruments. 

Data  Processinf: 

The  hourly  manual  data  is  averaged  over  the  month.  These  averaged  values 
are  plotted  versus  time  of  day.   7ero  points  representing  sunrise  and 
sunset  are  added  and  a  cuj^ve  is  drawn  tlirough  all  these  points.  This 
curve  thus  represents  a  monthly  average  insolation  curve.  The  curve  is 
integrated  numerically  to  yield  a  value  for  average  daily  insolation. 

Sample  curves  are  shown  in  Figure  7.  These  illustrate  the  interpolation 
procedure  for  stations  having  a  variety  of  data.  The  top  curve  is 
clearly  more  accurate  due  to  the  larger  data  base. 

The  average  solar  insolation  data  compiled  from  the  manual  stations  during 
September  and  October  are  given  in  Tabic  1.  Data  from  a  few  stations  was 
clearly  suspect  due  to  improper  orientation  of  the  transducer  or  reading 
error  and  is  omitted. 

Conclusion 

The  manual  portion  of  the  program  has  involved  a  large  number  of  students 
and   teachers  and  is  clearly  a  valuable  solar  educational  exercise.  The 
significance  of  the  data  will  be  evaluated  at  a  later  time  as  more  data 


7 

becQiics  available.      It  vjill  be  nccencary  to  compare  the  manual  d^i*, 
nearby '^rccordin,'^  ctationa  and  Kith  solar  noon  maximiimB  to  detem^iM 


statistically  neanincful  value  for  percentage  accuracy  or  reliabilftyT 
This  will  be  done  <-it   the  end  of  the  r.c}iool  year.     For  the  time  btfing,  ' 
the  manual  data  should  be  clasced  as  "pi'climinary  data"  because  the 
absolute  accuracy  has  not  been  detemiined.  < 
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Figure  7:     Examplec  of  I'onthly  Avcracc  Solar  Insolation 
from  Manual  Inatrument  Stations 


Avcraco  Daily  r>olar-Insolation 
/ . . .  1 Ly-\ 


*' 


Location                            „ 

\fZ  -nrv^i 

'm     /                                               '"SWaWMBBMI 

« 

ncpt 

Oct 

Anaconda 

1.94* 

3.03*             '^91 

Big  Timber 

Ko  Data  Available                     '  °^||H| 

Broadus 

5.31* 

5.10                ^W 

Chester 

3.89*                   >r,. 

Chinook 

4.24** 

5.83**                   ^''1 

Choteau 

4.63 

4.74 

Circle 

3.87* 

5.17* 

Columbus 

4.8? 

Deer  Lodge 

4.10 

5.15 

Ekalaka 

5.05** 

6.15 

Ennis 

?.83 

2.63* 

Port  Benton 

4.35** 

Hardin 

Ko  Data  Available 

Hobr.on 

/!.96** 

4.38** 

Jordan 

No  Data  Available 

1 

•  Lame  Deer 

Ko  Data  Av 

ailable 

Livincston 

4.91 

5.18      . 

Malta 

5.61 

Plentjn-food 

5.41 

Red  Lodce 

5.24** 

Roxindup 

4.11* 

5.39* 

Sidney 

3.03 

5.62* 

Superior 

3.04** 

Thompson  Palls 

4.70«* 

5.46** 

Tovmsend 

6.67** 

5.44** 

Vlhite  Sulphur  Spring! 

3 

4.02 

VJinnett 

5.87** 

5. 64** 

V?olf  Point 

5.05** 

5.79** 

»     only  one  daily  reading 
«*     more  than  4  daily  readings 


"*"     This  data  should  be  considered  preliminary  data  until   its 
accuracy  can  be  fully  assessed. 

Table    1:     Average  Daily  Solar  Insolation  Data  for  Manual 
Stationsj  September  and  October,    1977 
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APPIUDIX  Is 
TECHNICAL  DISCUSSION  OF  SOLAR  INSTRUT.IKNTATION  AND  ACCURACY 

Commorcial  Instruments 

Reviews  of  currently  used  Bolarimetors  and  pyranometers  are  given  by 
Blanco  (l),  Daniels  ( 2) ,  Haywood  ( 3) j  and  Hajumdar  (4).  Instruments 
based  on  thermo-electric  measurement,  calorimetric  measurement,  and 
photovoltaic  measurement  arc  available.  The  best  known  and  most  v?idely 
used  instruments  are  the  Eppley  and  Kipp  pyranometers.  The  Lppley, 
which  is  manufactured  in  the  United  States  and  commonly  used  here, 
consists  of  a  blackened  metallic  receiver  3  cm  in  diameter  surrounded 
by  a  white  or  silver  outer  ring,  A  copper-constantan  thermocouple  is 
attached  to  the  blackened  receiver  with  the  cold  junction  attached  to 
the  reflective  outer  ring.  The  temperature  difference  beiv/een  the 
black  body  and  the  reflective  body  gives  a  measure  of  the  radiation. 
The  reflective  ring  remains  at  nearly  ambient  temperature.  The  daily 
output  of  the  thermocouple  is  measured  by  a  recording  pyrcinometer.  This 
curve  ia  then  integrated  to  give  total  daily  radiation. 

The  receiving  surface  is  normally  placed  horizontally  and  covered  by  a 
hemispherical,  transparent  glass  or  plastic  cover.  The  direct  radiation 
as  well. as  the  diffuse  radiation  fall  onto  the  surface  and  are  recorded. 
Properly  calibrated,  the  Kppley  and  Kipp  pyranometers  have  an  accuracy 

Photovoltaic  Instruments 

Experiments  have  been  performed  to  test  the  applicability  of  silicone 
photovoltaic  cells  for  solar  radiation  measurements  by  Yellott  (5), 
Schoffer,  et.  al.  (6),  Banner,  et.  al.  (?),  and  l.'hillier  (8).  A  commercial 
instrument,  the  Sol-A-?^eter,  is  currently  manufactured  by  Matrix,  Inc.,  of 
Phoenix  under  a  license  from  John  Yellott. 

The  silicon  cell  consists  of  a  p-n   junction  of  ultra-pure  silicon.  Photons 
striking  the  cell  cause  electrons  to  be  displaced,  which  is  a  current.  Under 
proper  loading  conditions,  the  current  flow  is  directly  proportional  to  the 
incident  radiation.  This  behavior  makes  the  silicon  cell  a  candidate  for 

Bolar  measurement. 

The  advantages  of  the  silicon  transducer  are  its  relatively  high  electrical 
output  and  its  long-term  stability  and  reliability.  Its  disadvantages 
include  a  slight  temperature  dependence  and  a  spectral  sensitivity  which 
departs  from  black  body  behavior.  Extensive  experiments  are  reported  in 
the  references  above  which  compared  standard  Eppley  and  Kipp  pyranometers 
with  silicon  cell  pyranometers  (5f7»8).  These  experiments  show  that  the 
silicon  instruments  can  produce  results  within  t%.     Most  of  the  inaccuracy 
occurs  when  the  sun  is  at  low  angles.  The  variation  of  calibration  with 
optical  path  length  is  related  to  the  spectral  response  of  the  cell.  This 
effect  may  be  seen  in  the  calibration  data  for  the  SIMM  instruments  used 
in  this  program. 
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Instnunents  for  the  SI?'!!^  Pro/^am 

Instruments  based  on  the  silicon  cell  "transducer  were  chosen  for  use  in 
the  SIMM  program  for  two  reasonss      (l)   the  overall  cost  of  the  instruments 
was  minimized  and  (2)    the  characteristica  of  the  transducer  have  been  well 
documented.     That  is,   the  silicon  cell  measurements  are  related  to  the 
standard  radiation  instruraentB, 

This  second  factor  is  very  importajit.     There  are  a  large  number  of 
photovoltaic  and  photoresistive  transducers  available  at  the  present  time. 
Any  of  these  devices,  when  connected  to  the  proper  circuit  or  meter,  will 
give  Bome  change  in  output  with  variations  in  "light"  intensity.     The 
common  Cd-6  photographic  light  meter  is  an  example.     Due  to  linearity 
problems  and  spectral  response,   this  system  would  not  be  useful  for  energy 
measurements. 

The  immediate  purpose  of  the  SBIM  program  is  to  make  measurements  which  will 
accelerate  the  utilization  of  solar  energy.     It  is  therefore  important  that 
there  be  a  direct  relationship  of  the  measurements  made  to  the  measurements 
needed  for  solar  design.     Since  this  relationship  has  been  documented  in  the 
literature  (Ss^sTjS)?   no  additional  development  or  research  programs  are 
needed  to  establish  "the  meaning"  of  the  measurements. 

Prof.  John  Yellott  has  U.S.    patents  covering  some  aspects  of  solar 
measuring  instruments  using  silicon  cells.     I  have  described  the  SUM 
progrsun  and  the  instruments  to  him.     He  has  assured  me  that  the 
Sim-I  instruments  (l)   would  not   infringe  on  his  patents  and  (2)  would  be 
quite  accurate  and  easy  to  calibrate. 

Description  of  the  SBUI  Instrument  Package 

The  instrument  package  distributed  to  the  high  school  measuring* stations 
consists  of  two  solar  radiation  measuring  units,    the  SIKK-1    and  the  SILM-IR. 
Photographs  of  the  instruments  are  shovm  in  Figure  1-1, 

SIM-1    is  a  hand-held  instmiment.      It  consists  of  a  silicon  cell  mounted  on 
a  tiltable  plate.     The  plate  is  hinged  to  a  small  instrumentation  box  which 
contains  the  indicating  meter.     The  instrument  is  self^ovjercd  and  needs  no 
batteries.     It  is  used  in  this  program  as  a  back-^ip  instrument  and  as  an 
educational  tool. 

The  Sn-m-IR   instrument   incorporates  a  strip  chart  recorder  so  that  a 
continuous  record  of  solar  energy  can  be  obtained.     The  recorder  drive  is 
powered  by  100  vac.     The  transducer  is  separate  from  the  recorder  and 
consists  of  a  silicon  cell  mounted  on  an  angled  plate.     This  assembly  is 
designed  for  attachment  to  an  exterior  wall  or  roof.     The  recorder  is 
located  inside  the  building  on  a  convenient  shelf  or  cabinet. 

Tramsducer  Orientation  of  SII-TI1-1R 

The  silicon  cell  in  the  SBM-IR   transducer  will  be  mounted  facing  south 
and  at  an  angle  of  60  degrees  from  the  horizontal.     This  orientation  is 
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chosen  because  it  is  nearly  the  optimum  winter  collection  angle  for 
Montana^  ,^The  data  taken  at  this  angle  can  accurately  be  transformed  to 
other  similar  angles*  V/eather  bureau  data  recordss  radiation  on  a  horizontal 
eurfacco  Transforming  this  data  to  practical  collector  angles  produces  error, 
(in  Montana  in  the  winter  the  altitude  of  the  sun  is  about  20  degrees.)  It 
is  the  philosophy  of  this  program  to  make  measurements  which  will  implement 
the  utilization  of  solar  eneror.  This  goal  is  achieved  by  taking  the  data 
at  orientation  matching  the  ultimate  device. 

Calibration  and  Error 

Calibrations  of  the  instruments  were  done  by  comparing  the  output  of  the 
prototype  SBM  instrument  with  an  Eppley  pyi-anomcter  and  a  Kipp-^oncn   i;  ■ 
pyranometer.  These  calibrations  were  performed  with  the  SIH-I  inctrument 
in  the  horizontal  position,  to  match  the  orientation  of  the  Kpploy  and 
Kipp,  Eppley  and  Kipp  instruments  are  probably  the  most  common  instruments 
used  for  solar  enerfsf  measurement  around  the  world  and  their  performance 
is  well  understood,  (See  Figure  1-2,) 

The  results  of  a  calibration  run  for  an  entire  clear  day  are  shown  in 
Figure  1-3.  The  calibration  results  show  that  the  SIMM  can  be  calibrated 
to  within  a  deviation  of  IJj  during  the  middle  part  of  the  day.  In  the  early 
morning  and  the  late  evening  the  deviation  rises  to  5-'\%»     This  phenomenon 
is  explained  by  the  spectral  response  of  the  silicon  cell  and  is  well 
documented  in  the  solar  literature.  These  errors  at  the  beginning  and  end 
of  the  day  are  not  as  large  as  they  may  seem  because  the  absolute  energjr  is 
low* 

The  error  in  total  energy  on  this  particular  day  is  shovm  in  Figure  11  to 
be  less  than  l^^j.  Because  the  absolute  accuracy  of  the  Eppley  and  Kipp 
instrrnnents  is  only  about  %f   it  can  be  seen  that  the  SIMM  instrument  is 
well  within  this  band  of  error. 

Errors  in  Data  Reduction 

Data  is  read  from  the  strip  charts  by  hand  by  a  Single  staff  member.  The 
average  solar  radiation  over  a  IS-minute  time  interval  is  recorded  as  a 
two-Hiigit  number.  Approximately  60  such  numbers  are  read  from  each  chart. 

Some  technical  judgment  is  required  to  read  the  charts  on  days  when  there 
is  mixed  overcast.  These  curves  are  jagged  and  the  rough  curve  must  be 
averaged  over  the  time  interval. 

To  test  for  reading  variations,  10  charts  were  selected  at  random  from  our 
files.  These  charts  were  re-processedj  by  the  same  staff  member,  and  the 
results  compared  to  the  results  from  the  first  readings.  The  results  of 
this  test  are  shown  in  Figure  1-4  It  can  be  seen  that  the  average 
uncertainty  in  the  data  reading  process  is  2,45^ 

Errors  in  Data  Handling 

These  numbers  (the  raw  data)  are  put  into  the  computer  system  through  a   ' 
remote  terminal.  As  the  numbers  are  put  into  the  system  the  computer 
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checks  for  olrviouc   errors— larf^e  numborr.j   ncf^itivc  niun"bersj    etc.     After 
the  numbers  from  a  maximum  of  14  days'   charts  are  entered,   the  computer 
is  placed  in  a  LIST  mode«     The  numbers  are  printed  out  and  the  operator 
checks  them  a£:ainGt  the  original  numbers^     Errors  in  punching  are  then 
corrected  in  an  EDIT  mode.     The  final,   corrected  data  is  then  stored  in 
the  computer*     The  possibility  of  sifTiificant  error. in  this  stage  of  the 
data  manipulation  is  negligible,  due  to  the  two-fold  error  checking 
procedure. 
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Figure  1-4:  Data  Reduction  Reliability 
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Fig-are  1-3: 
Calibration   of  SllZ'i   Instruments 
Results  of  a  Calibration  Run   for  an  Entire  Clear  Da;y 
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APP'IfDIX  II: 
DATA  niiDUCTIOM  PitOCICBUlilS 

Introduction 

The  general  flovf  chart  for  the  data  reduction  is  nhown  in  Pi{jureII-1. 

The  solar  energy  data  concicts  of  a  continuouD  record  of  solar  radiation 
versue  time.     At  any  instant  of  time  the  amplitude  of  the  chart  can  bo 
multiplied  by  thCpCalibration  factor  of  the  inEtruinent  to  yield  the  solar 
radiationj  BTU/ft  --hr  (or  H' /m  ).     A  sample  set  of  curves  for  Helena  is 
shovm  in  Figure    II~?„ 

One  type  of  data  that  ic  uaeful  is  the  total  radiant  energy  falling  on 
the  trancducer  durin^j  the  entire  day,   yL'-iirr./m  -day.     Tliis  value  can  be 
determined  from  the  area  under  the  chart  record.     The  appi^oximate 
integration  of  the  experimental   curve  is  performed  by  constructing;  an 
approximate  curve  ujnd  then  intcp:'atin2  that  curve.     The  approximate 
curve  consists  of  values  averaj^jed  over  an  inter/al  of  15  minutes.     The 
integral  of  the  approximate  curve  is  then  simply  the  sum  of  these 
incremental  average  values  through  the  day. 

The  accuracy  of  this  method  of  integration  depends  on  the  determination 
of  the  average  values.     The  accuracy  will  vary  from  record  to  record, 
depending  on  the  exact  character  of  the  rr'.diation,     A  Gmootli  cui-vo,   cuch 
as  from  a  clear  da^,",   can  be  more  accurately  integrated  than  the  irregular 
curve  from  a  cloudy  day.     Because  60  or  more  increments  are  used  to 
approximate  each  curve,    errors  in   judgment  v.'ill  tend  to  cojicel.     A 
specific  evaluation  of  the  reliability  of  this  procedure  is  given  in  a 
later  section. 

Comrtuter  CutT)ut 

Figure  II-3  depiits  tlBdata  reduction  process  for  one  day.     The  basic 
solar  insolation  curve  is  shown  in  Fi^^iire   3a.     The  set  of  numbers 
representing  the  insolation  curve  is  shown   in  Figure   3b.     This  set  of 
numbers  is  chocked  for  accuracy  and  then  used  to  calculate  the  total 
and  threshold  energiecj  Figure   3c» 

Results 

Dxanples  of  the  final  daily  radiation  data  for  one  station  (Billings) 
for  12  days  of  operation  is  shovm  in  Figure  II -4 (Days  are  numbered 
consecutively J  starting  from  January  1  as  Day  1.)  This  data  is 
averaged  to  obtain  the  values  Bhov;n  in  the  final  compilation. 
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Figure  IW  :  plow  Chart  for  Solar  Data 
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irUI  ACTlVITIli;  DUiaNG  FlllM,  OU..uiT"H 


Avards 

The  volunteer  teachers  vere  unifomly  cooperative  and  enthuciastic 
during  the  first  year's  operation  of  Hlini,  They  were  generally  prompt 
in  sending  in  the  solar  data  and  helped  on  come  occasiono  to  make 
minor  adjustments  in  the  recorder  and  transducer. 

To  reco^Tiize  their  contribution  we  made  up  "Certificates  of 
•Appreciation."  The  authority  of  these  certificates  wao  amplified 
with  the  signatures  of  John  Orth  and  Govcx'nor  Judjrrc,  A  sample 
certificate  is  shovm  in  Figure  5» 

Manual  Data  Stations 

V.'e  have  "been  receiving  data  from  the  30  manual  stations  since  the 
beginning  of  the  school  year.  L'e  requested  that  the  utudcntc  or 
teachers  take  periodic  solar  intensity  mcaf-aremcnts  usin^  the  manual 
instruments  and  enter  them  on  a  chart.  Tlicsc  charts  are'^sent  to  us 
monthly. 

Our  results  have  been  variable  and,  ac  expected,  depend  on  the 
dedication  of  the  teachers  and  students.  V.ozt   of  the  stations  arc 
returning  daily  noon  readings  as  v.-e  requested.  About  25-30;'  of  the 
stations  go  far  beyond  our  remiest  and  send  us  hourly  data.' 

Ve  are  just  beginning  to  compile  this  data.  It  v/ill  be  necessary  to 
do  a  careful  statistical  analysis  of  this  data  before  we  are  prepared 
to  include  it  v,'ith  our  data  from  the  automatic  recordinfi;  stations. 

I  am  convinced  that  much  of  the  manual  data  v;ill  have  cicnificjijice 
when  properly  interpreted.  The  vol-LUitecrs  are  certainly  accraing  a 
lot  of  educational  value  from  their  participation.  The  data  from  the 
manual  stations  is  presented  and  discussed  in  the  data  section  of  this 
report. 

Direct  Data  Reduction 

Our  present  data  reduction  process  (Appendices  I  and  II)  involves 
considerable  handling  of  tlie  data.  The  numbers  arc  read  from  the 
chart  J  written  dovm,  entered  into  the  computer,  and  checked.  IJach  of 
these  steps  must  be  done  for  an  average  of  over  1000  numbers  each  day. 
This  is  time  consuming  and  the  multiple  steps  expose  the  data  to 
possible  errors. 

Initial  v;ork  has  begun  on  a  curve  reader.  A  prototype  is  shorn  in  the 
photograph  in  Figure  II t-1 .  The  purpose  of  this  device  is  to  streamline 
the  data  reduction  process.  The  operation  consists  of  the  follov;ing 
steps J 
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(1)  The  chart  paper  containing  the  rolar  insolation  curve  sent 
by  the  high  school   in  placed  in  the  chart  drive  mechanism* 

(2)  The  calibration  factor  for  the  particular  recorder  ic  net 
on  the  controlo  on  the  ri^jht  h.-md  part  of  the  inctrument. 

(3)  The  r.ero  of  the  chart   is  adjusted. 

(4)  The  chart  drive  is  activated  and  the  pointer  is  manually 
positioned  to  follow  tho  solar  insolation  cuivc. 

The  output  of  the  curve  reader  is  a  time-viirying  electrical  volta^je 
which  is  aji  analpc  of  the  solar  insolation.     The  time  base  of  the 
analog  signal  is  about  one  minute  and  represents  about  12  hours  real 
time.     Once  the  instrument  has  been   zeroed  the  entire  tv;o  v;eeks  of 
chart  are  read  sec[uentially. 

This  data  si^al  will  be   entered   directly  into  a  computer  for  integration, 
processinff,    and  data  storai^e.     This   system  vjill   reduce  the  amount  of  tine 
nov;  needed  to  read  the  curves  manually.     There  v;ill  also  be  less 
opportunity  for  accidental    error  since  the  data  sifpnal  will  be  entered 
directly  instead  of  manually  into  tlic  computor„ 
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APPniDIX  TV: 
COlJCULTAin'S 

Listed  below  aro  consultants  \^lao  have  contributed  to  the  SIITIl  Procraxn. 
A  brief  statement  ic  included  to  dcncribc  their  areas  of  expertico. 

John  Yellott 

K  B.S.  I  O.S.J.  I  Registered  Mechanical  Lncincer  in  Arizona,  Ilai^'land,   and 
New  Jersey?  Fellow,  A.S.K.P..,   A.S.H.K-A.E.,  A.A.A.S.  Ari7-ona  AcadoRiy  of 

Science.  r.;^;,  V  ' 

Yellott  SolarllnerQr  Laboratoiy,  901  V?est  ?.l  Caminito,  Phoenix,  Arizona 
65021. 

Ur.  Yellott  is  a  recognised  expert  in  the  areas  of  colar  incolaiion 
measurement  and  nolar  enercy  utilisation.     His  profcsoional  activities 
mi  publications  would  rectiiire  several  pa^cs  of  description  and  are  thus 
not  included  here.     lie  has  published  orifjinal  papers  and  holds  patents 
in  the  area  of  solar  insolation  measurement  using  silicon  cell 
transducers,   the  type  of  device  used  in  tlie  5II3i  Procram. 

Mr.  Yellott  advised  us  durin^r  the  early  staces  of  the  SIID'i  Procran 
©^nceming  instra-ncnt  development,     lie  has  reviewed  our  calibration 
procedures,   instrument  dcsicn,   and  data  processing.     His  advice  has  civen 
us  valuable  assurance  of  the  level  of  accuracy  of  our  measurements.     Ve 
feel  fortunate  that  an  internationally  known  solar  expert  is  associated 
with  the  sum  Program. 

Bob  Blair 

15e  iS  o  Li » y^  9  ' 

Department  of  Aeronautics,  University  of  V'ashinctonj  Seattle,  V.'ashinerton. 

llr.  Blair  is  en   electrical  enfp.ncer.  He  has  had  broad  experience  in 
desicninc  an<i  building  instrumentation  systems  for  research  and  testing 
applications,  including  the  Boeing  Vind  Tunnel,  KED  and  plasma  research, 
shock  tubes,  holography,  experimental  mechanics,  etc. 

Mr.  Blair  helped  to  locate  reliable  and  economical  recording  instruments 
for  the  Slim  Program.  He  continues  to  give  us  expert  opinions  on 
electronics  questions  regarding  data  reduction. 

Phillip  Mast 

Ph.D. 

Scientist,  Naval  Research  Laboratory,  V'ashington,  B.C. 

Dr.  Hast  has  had  over  20  years  of  research  experience  in  the  areas  of 
materials,   instnamcntationg   and  computer  applications.     His  experience^ 


ranees  from  ori^;inal   contributions  to  the  thooricc  of  infoimation  to 
practical  desi^Ti  and  construction  of  computer  liardv;are. 

Dr.  I'aot  has  advised  us  concerning  our  data  procosfiin,'^  methods.  Our  f;oal 
is  to  incorpoi'atc  cuid  utilize  statc-of-thc-art  computer-assisted  routines 
for  proccBsin;^  the  jolar  insolation  data.  By  ucinc  these  methods  v/e  hope 
to  maximize  the  accuracy  and  minimi :;c  the  cost  of  our  data  processing;  and 
storaso, 

Leirry  Thomas 

Science  and  Mathematics  Supervisor,  Office  of  the  Superintendent  of  Public 
.Vjlnstruction,  Helena,  Ilontana     59^0'' • 


I^Currcntly  teaching  in  Helena.) 

it. 

Mr.  Thomas  contributed  a  lot  of  time  and  ideas  during  the  early  Gt&.,[:es  of 
rthe  SIi>!  Program.  Ko  helped  to  contact  the  }iic:h  school  volunteers  arid 
:,wrotc  many  letters  explaining  ^■'ho  ve  were  and  what  we  planned  to  do»  His 

advice  and  help  in  the  communications  area  vjas  invaluable. 

Mr.  Thomas  consulted  v;ith  us  concemin-j  the  content  of  the  "nandboo'<:  of 
Solar  Lnergy  Experiments,"  Pie  read  the  first  edition  of  this  handbook  and, 
as  an  expert  in  science  teachinr,  pave  his  advice  on  changes  in  format  and 
content  of  the  handbook,  llis  suggestions  vdll  be  incorporated  in  future 
editions. 
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THE  SPACING  OF  SOL/J?  MEASUREMENT  STATIONS 


ScieniiBts,   climaiologists  and  ma-themeticians  have  worked  on  the 
question  of  maximum  allowable  inter-sta-tion  spacing.     The  hasic  goal 
is  to  distribute  solar  insolation  recording  stations  close  enou^ 
together  bo  that  values  of  solar  insolation  between  the  measuring 
Btations  can  be  estimated  -with  a  mininrora  or  at  least  acceptably  email 
error.     This  .error  is  called  "interpolation  error"« 

Some  of  the  most  complete  work  in  this  area  is  due  to  L..-S.  Gandin 
ia  The  Planning;  of  Meteorological  Station  Networks  from  the  World 
Meteorological  Organization  of  the  United  Nations.     Other  important 
material  is  available  from  D.  0.   E,   contract  EG-77-S-<> 5-5362,  Solar 
Radiation  Data  Sourcest  Applications;,  gjid  Network  DesiE?i»     Statistical 
argoments  and  data  from  the  latter  are  reproduced  in  detail  in 
Appendix  I, 

To  summarize  J  Gandin  analyzed  48  stations  and  by  a  statistical  argument 
concluded  that  for  an  optimal  interpolation  error  of  0.11   the  inter- 
station  distance  should  be  about  15O  kilometers  (91  miles).     Figure  la 
shows  such  an  optimal  grid  covering  the  state  of  Montana.     This  grid 
contains  a  total  of  26  stationso     For  comparisonj  Eigure  lb  shows  the 
actual  distribution  of  current  SDttl  measuring  stations. 


Conclusions 

The  current  30  station  SIMM  network  is  very  similar  to  Gandin's 
"optimal"  26  station  network.     Since,   at  any  given  time  there  may  be 
1  to  3  SH'M  stations  off  line  due  to  operator  error  or  machine  mal- 
funcition,   the  current  30  station  network  is  surprisingly  close  to 
Gandin' B  "optimal"  network.     Having  a  few  extra  stations  allows  the 
SBM  network  to  be  tolerant  of  these  errors  and  thereby  has  a  built-in 
safety  factor. 

It  appears  that  the  current  30  station  grid  has  a  satisfactory  density  at 
present  and  should  be  maintained.     Expansion  of  the  grid  would  possibly 
be  justified  in  some  mountain-valley  regions  of  the  state  where  there 
may  be  significant   local  weather  variations  (Bozeman  and  Belgrade,    for 
example).     The  current  network  has  its  highest  density  in  the  western 
and  central  region  and  this  appears  to  be  justified. 


(a)  'theoretical  Optiffltira"26  Station  Grid 


O 
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(b)  Current  30  Station  Grid 
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HOW  MANT  YEARS  OF  SOLAR  DATA  ARE  NEEDED? 


A  basic  concept  of  data  collection  Ib  illuBtrated  in  the  conliaence 
curve  Bnovm  in  Figure  2,  While  literally  thousands  of  pages  in  technical 
hooks  and  "papers  address  this  basic  pxxthlem  in  a  detailed  coverage, 
this  material  is  beyond  the  scope  of  this  proposal.  In  this  area  we 
rely  primarily  on  the  expert  judgment  of  scientists  working  in  the 
field  of  solar  radiation  measurement  rather  than  deriving  the  statietical 
arguments  in  detail. 

For  purposes  of  discussion  it  is  helpful  to  consider  the  ^'confidence 
curve"  in  Figure  2.  This  curve  illustrates  that  when  no  data  is 
available  the  confidence  level  is  "0".  The  first  data  points  result 
in  a  large  increase  in  confidence  above  the  "0"  level  (the  first  year 
SIMM  measurements) «  The  later  part  of  the  curve  tapers  off  illustrating 
the  smaller  increases  in  confidence  as  larger  amounts  (greater  thaa 
10  years)  of  data  points  are  available. 

There  are  two  important  reference  points  on  this  curves  at  two  years 
and  ten  years.  The  "two  yeaj"  point  represents  the  position  of  the 
second  year  SlMi   programf  that  isj  at  the  end  of  this  current  grant 
period.  At  the  end  of  the  two  year  milestone^  the  Department  of  Energy 
will  consider  including  SIM   data  in  the  national  solar'  energy  data  bank. 
They  feel  that  solar  data  over  time  spans  less  than  two  years  is  not 
significant. 

The  other  milestone  is  the  10  year  i>oint.  The  NOAA  National  Wsather 
Service  requirement  for  publishing  average  solar  energy  data  is  10  years 
of  continuous  solar  measurements.  Many  meteorologists,  incidentally, 
favor  an  11  year  minimum  period  for  establishing  average  solar  energy 
data*  This  preference  results  from  the  fact  that  the  natural  solar 
cycle  is  11  years « 

These  reference  points  at  2  years  and  10  years  shed  some  li^t  on  the 
(luestion  of  accepted  time  intervals  for  securing  representative  solar 
data  and  help  to  evaluate  the  status  of  our  current  measurment  program. 

Additional  informa-tion  results  by  reviewing  some  actual  SIMM  data  in 
the  next  sectiona 


Annual  Data  Flicjaatlone 

There  exist  naturgil  variations  in  solar  energy  from  year  to  year.  These 
fluctuations  are  msdnly  due  to  fluctuations  in  local  cloud  cover.  The 
weather  pattern  each  year  is  to  some  extent  unique.  If  data  is  taken 
for  a  long  period  to  time,  these  yearly  fluctuations  can  be  distinguished 
from  overall  predominant  weather  patterns.  It  is  these  overall 
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characteristics  that  are  of  practical  use  in  solar  energy  ntilizatione 

A  solar  heating  systemf   for  example^   is  expected  to  operate  over  the 
lifetime  of  the  Iruilding  which  may  be  20  to  5°  years.     The  solar  Bystem 
muBtj   thereforej  he  Eized  to  correspond  to  some  ''average"  radiation 
existing  throughout  this  20  to   50  year  time  period.     A  rule  of  thumh 
for  solar  desi^  aBsumes  that  annual  fluctuations  in  solar  radiation 
will  he  normally  contained  in  a  hand  of  t  10-1^  around  the  long  time 
(10  years  or  greater)   average* 

Solar  data  fluctuations  for  the  first  four  months  of  the  years  1977  and 
1978  for  14  statione  has  heen  plotted  in  Figure  3.     These  deviations 
are  calculated  hy  tsQ;:ing  the  difference  in  year  1977  and  year  1978  data. 
The  average  daily  total  radiation  for  each  month  was  used  as  the  hasis 
of  the  graph.     A     ±  10^  and    ±   15^  hand  are  also  plotted  on  the  graph 
to  compare  "expected"'  fluctuations  to  measured  fluctuationse 

Examination  of  the  gi-aph  shows  that  there  are  raore  points  outside  the 
t  1C^  hand  than  there  are  within  the  hand!     This  evidence  indicates  that 
at  least  one  of  the  years  is  prohahly  not  a  "typical"  or  "average"  year. 
A  second  observation  of  Figure  3  ib  that  if  the    ±  105J  fluctuation  hand 
were  displaced  upward,   more  points  would  he  contained  in  the  hand.     This 
could  indicate  that  (1)   the  solar  energy  during  1977  was  unusually  high 
or  (2)   the  energy  during  I978  was  unusually  low.     Sufficient  data  is  not 
available  to  scientifically  or  precisely  answer  this  question  at  this 
time.     Some  reference  can  be  made  by  looking  at  NOAA  data. 

Comparison  with  Great  Falls  and  Glas^w  NCAA  data 

Long  term  averages  (IO  years)   of  solar  energy  falling  on  a  horizontal 
surface  are  available  from  tables  compiled  by  the  National  Weather 
Service  of  NOAA  for  (only)   two  locations  in  Montana;     Great  Falls  and 
Glasgow.     SIMM  measurement  stations  were  located  at  these  two  cities 
for  the  main  purpose  of  providing  a  comparison  of  our  current  short-tTsnn 
to  available  long-terra  solar  data.  The  NOAA  average  horizontal  data  was 
transfomed,  using  accepted  ASHRAE  methods,   to  a  60°  surface  for 
comparison  to  SIMM  data. 

Figure  4  shows  a  graph  containing  data  from  SIMM   1977  and  SIMM  I978 
along  with  long-term  average  NOAA  data  for  Great  Falls  and  Glasgow  for 
the  first  six  months  of  the  year.     There  are  several  conclusions  that 
can  be  drawn  from  this  comparisons 

0       For  1977  and  1978  the  solar  energy  available  in  Great 

Palls  and  Glasgow  during  the  months  of  Januaiy-j  February 
and  March  averaged  40^  lower  than  the  10  year  average. 

(D        The  solar  energy  in  1978  for  Great  Falls  and  Glasgow  was 
below  average  for  the  entire  spring  heating  season. 
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Figure    3  :  Yearly  deviations  in  monthly  average, 
daily  solar  energy  for  20  stations 
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d)  The  Bolar  energy  for  1977  i^   Glasgow  and  Great  Falls 
during  the  months  of  April,  Hay  and  June  was  "normal*'! 
i.e.  it  fell  within  the  t     1C^  variance  "band  when 
compared  to  long-term  average  values* 

Conclusion 

Solar  data  spanning  a  longer  time  period  is  needed  to  "build  a  firm 
solar  energy  'data  base  for  Montana.  The  evidence  Indicates  that 
"both  1977  and  1978  have  had  an  unusually  small  amount  of  solar  energf 
during  the  crucial  heating  season  in  mid-winter  and  are  not  typical* 
Solar  utilization  systems  designed  to  this  short-term  data  would  tend 
to  "be  ±00  large  and  hence  less  economically  feasi'ble. 

The  confidence  curve  in  Figure  2  shows  that  a"bout  10  years  of  data 
is  needed  to  provide  a  90^  confidence  level.  This  would  "bring  the 
average  solar  energy  values  approximately  within  the  ±  10?^  fluctuation 
"band.  After  this  10  year  point  the  pay-off  in  accumulating  more  data 
is  small.  For  times  less  than  10  years  each  additional  year  of  data 
yields  a  significant  increase  in  confidence.  It  appears  from  the  curve 
that  a  practical  minimum  time  span  for  this  measurement  prograaa  would 
"be  5  or^  6  years. 

On  August  4g  1978  I  was  invited  to  Oregon  State  University  to  speak  on 
our  SIM^i  solar  measurement  program.  The  program  was  attended  lay 
meteorologists  and  scientists  from  5  northwest  states  who  are  currently 
measuring  solar  energy  resources  and  are  planning  expansions  of  these 
programs.  Diiring  my  talk  I  asked  for  opinions  from  this  group  on  the 
length  of  time  required  to  establish  valid  "solar  energr  averages". 
The  general  opinion  was  that  5  to  10  years  of  continuous  measurement 
would  be  a  practical  goal. 
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SOLAR  ENERCn  BATA  FROM  SBM  PROGRAM 


Introduction 

This  section  descrjloes  the  data  handling  in  the  SDIM  progrsuno     The 
present  status  of  the  data  for  1973  is  indicated  "try  the  bar  graph  in 

Figure  5  » 

pata  Reduction 

The  general  flow  chart  for  the  data  reduction  is  shown  in  Figure  6  » 


The  solair  enev^   data  consists  of  a  continuous  record  of  solar 
radiation  versus  time.  At  any  instant  of  time  the  amplitude  of  the 
chart  can  he  multiplied  hy  the  calibration  factor  of  the  instrument 
to  yield  the  solar  radiation^  Btu/ft^-  hr  (or  ln;/m^  )»  A  sample  set 
of  curves  for  Helena  is  shown  in  Figure  7  «  Figure  8  shows  curves 
from  four  different  cities  taken  on  the  same  day* 


One  type  of  data  that  is  useful  is  the  total  radiant  energy  falling 
on  the  transducer  during  the  entire  day,  Idf-hrs/m'^  -  day»  This  value 
can  be  determined  from  the  area  under  the  chart  record. 


Reading  Cheirt  Records 

The  solar  energy  data  is  recorded  and  sent  to  our  office  in  the 
form  of  strip  chart  records*  Last  year  we  read  numbers  from  these 
charts  and  then  enetered  them^  by  hand,  into  a  computer  system." 
Although  this  method  worked  quite  well,  a  lot  of  time  was  recjuired 
to  process  the  data.  The  associated  computer  rental  time  was  also 
expensive. 

A  study  of  this  problem  indicated  that  methods  should  be  developed 

for  entering  the  chart  data  directly  into  the  computer  using  a 

curve  reader,  Devices  for  inputting  graphical  data  have  been  available 

for  several  years  but  their  cost  has  been  prohibitive  for  a  project 

of  our  scale  and  budget.  Cost  reductions  and  developments  dxiring  the 

past  few  months  in  the  computer  field  have  placed  usable  systems  within 

reachn. 

We  now  have  in  our  office  a  32  K,  8  bit  digital  computer  system  with 
dual  magnetic  disc  stora;je.  This  system  has  a  keyboard  and  video 
terminal  for  communication  with  the  computer.  A  special  feature  is 
a  gi'aphic  input  tablet  which  is  interfaced  directly  to  the  computer. 


c:vATU':  OF  1970  gol/j:  dv:'a  icun. 
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The  input  tablet  consists  of  a  flat  platen  about  12  inches  square  and 
a  cursor  which  looks  like  a  pen.     VJhen  the  pen  is  touched  to  the  platen, 
the  Bystem  outputs  the  i-y  coordinates  of  that  point  on  the  platen 
directly  into  the  computer  in  digital  form.     Thus  when  a  curve  is 
continuously  traced  on  the  platen,   a  series  of  x-y  coordinates  which 
represents  the  curve  is  siimiltaneously  input  to  the  computer. 

If  a  solar  energy  chart  record  is  placed  on  the  platen,  the  graphical 
information  can  be  directly  entered  into  the  computer  by  simply 
tracing  the  curve  on  the  chart  record.     This  device  eliminates  nmch 
of  the  "hand  work"  and  promises  to  improve  accuracy  as  veil  as  econocny, 
A  photograph  of  the  system  is  Ehoim  in  Figure  9, 

New  Data  Reduction  Process 

A  schematic  flow  chart  of  the  data  reduction  and  storage  using  this 
system  is  sho'rfn  in  Figure  10.     The  main  feature  of  this  process  are  the 
use  of  the  graphics  input  tablet  to  enter  the  solar  energy  chart  records. 
The  reduced  data  is  stored  on  magnetic  floppy  discs  in  the  form  of 
time-of-day  averages  of  solar  intensity.     The  file  structure  of  the 
data  storage  is  outlined  in  Figure  11, 

Once  the  solar  data  is  in  computer-readable  form  on  the  magnetic  disc 
it  can  be  easily  and  efficiently  processed  and  ouput  in  aiay  desired 
format.     A  primary  output  is  Tables  of  Average  Insolation  in  the 
format  used  during  1977  in  the  SBM  program.     Direct  computer  simulations 
of  solar  systems  woxild  also  have  been  run  using  this  system. 

Threshold  Radiation 

In  addition  to  total  radiation  data,  we  are  also  calculating  the 
radiation  above  a  series  of  thresholds.  The  reason  for  calculating 
threshold  data  derives  from  the  operating  characteristics  of  nearly  alJ. 
solar  collectors.     A  solar  collector  is  nonnally  operating  only  during 
periods  of  relatively  high  solar  intensity.     This  is  due  to  the  losses 
inherent  in  any  collector.     In  early  morning  or  evening  hours  or  during 
other  periods  of  low  radiation,  the  collector  will  lose  more  heat  than 
it  will  gain.     The  radiation  falling  on  the  collector  must  be  of  an 
intensity  great  enough  to  raise  the  collector  surface  temperature  ABCfVE 
the  temperature  of  the  storage  reservoir. 

For  example,   if  the  storage  reservoir  has  a  temperature  of  125°  F  and 
the  ambient  air  is  25°     and  if  the  collector  is  "double-glazed  ,  the 
heat  loss  through  the  glazing  will  be  about  0.68  Btu/^t'^  -  hr  -      F« 
or  about  68  Btu/Oir-ft^  . 

Assume  that  other  conduction  losses,  when  added  to  this,   result  in  a 
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to-tal  loss  of  90  Biu/hr-f-t2  ^       Thus,   if  tho  radiant  in-tensiiy  is  less 
than  this  amount,   the  collector  will  operate  at  a  net  LOSS. 

A  series  of  threeholds  from  0.1  kw/m^  to  1.2  w/m^  ,   spaced  0.1  KW/m^ 
is  xised.     Given  the  operating  threshold  of  a  particular  collector,   -it'is 
possible  to  interpolate  the  data  to  detemine  the  actual  energy  available. 

Data  Tables 

To  demonstrate  some  of  the  data  that  is  currently  on  record  we  have 
printed  the  data  tables  in  Appendix  H.   The  format  of  the  tables  is 
similar  to  the  format  used  for  the  tables  presented  in  last  year's 
final  report.     The  day  numbers  refer  to  the  day  of  the  year  with 
January  1  being  day  1,   etc.     Individual  daily  totals  as  well  as  a 
breakdown  of  the  intensity  spectrum  are  given. 

Fault  Codes 

As  we  enter  the  data  we  also  record  a  "fault  code»«  for  the  data.  If  the 
data  is  missing  or  suspect  due  to  instrument  or  operator  malfunction 
we  record  that  data  on  the  disc  by  entering  the  appropriate  fault  code 
number.  This  will  allow  us  to  scan  through  the  data  at  year's  end  and 
compile  and  categorize  the  frequency  and  types  of  data  problems  or  losses. 

Some  small  faults  in  the  data  are  recoverable  and  we  have  routines  built 
into  the  computer  program  to  adjust  these  curves  so  that  we  can  still  use 
them.  Any  data  that  we  have  adjusted  is  noted  vn.th  an  appropriate  fault 
code  so  that  we  can  separate  it  from  normal  valid  data.  In  the  tables 
recovered  data  is  noted  by  parentheses  around  the  total  energy  value  for 
that  day,  . 


Accuracy  Checks 


We  are  currently  xnmning  accuracy  checks  by  repeating  about' 300  data 
days  of  chart  records.  The  comparison  of  these  records  will  give  us  a 
good  idea  of  our  repeatability.  VJe  plan  to  do  periodic  random  checking 
of  the  data  entry  process  to  maintain  a  running  record  of  the  accxiracy 
of  our  data  reduction  process.  Current  test  results  show  that  our  average 
accuracy  is  "better  than  2.5^. 


Conclusion 


It  is  significant  to  note  that  our  data  reduction  system  closely  parallelt 
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an  elaborate  Bystem  devoloped  imder  a  D«  0»  E»  contract  for  reclaiming 
Bolar  energy  ctart  data  in  the  United  States,  We  can  thus  feel  confident 
that  we  are  technically  "up  to  speed"  in  this  area  of  reducing  Eolar 
energy  charts.  An  economic  comparison  in  another  section  of  this 
proposal  shows  that  we  are  doing  this  data  reduction  jol)  for  ahout  20^ 
of  D.  0.  £.   costs. 
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TECHNICAL  DISCUSSION  OP  SOLAR  INSTRTOffirJTATION  AND  ACCURACY 

Cormnercial  Instniments 

Reviews  of  currently  used  solariraeters  and  pyranometers  are  given  "by 
Blanco  (l),  Daniels  (2),  Hayvood  (3)i   and  MajTimdaj'  (4)*     Inst  ram  ents 
"based  on  thermo-electric  measurement j    calorimetric  measurement j   and 
photovoltaic  mcaBureiiient  are  availahle.     The  "best  known  and  most  widely 
used  instilments  are  the  Eppley  and  Kipp  pyranometers.     The  Eppleyj 
which  is  manufactured  in  the  United  States  and  commonly  xised  here^ 
consists  of  a  blackened  metallic  receiver  3  cm  in  diameter  surrounded 
"by  a  white  or  silver  outer  ringo     A  copper-  constantaa  thermocouple  is 
attached  to  the  "blackened  receiver  with  the  cold  junction  attached  to 
the  reflective  outer  ring.     The  temperature  difference  "betvreen  the 
"black  "body  and  the  reflective  "body  gives  a  measure  of  the  radiation. 
The  reflective  ring  remains  at  nearly  am'bieni  temperature.     The  daily 
output  of  the  thermocouple  is  measured  "by  a  recording  pyranometer.     This 
curve  is  then  integrated  to  give  total  daily  radiation. 

The  receiving  surface  is  normally  placed  horizontally  and  covered  "by  a 
hemispherical,  transparent  glass  or  plastic  cover.     The  direct  radiation 
as  well  as  the  diffuse  radiation  fall  onto  the  surface  and  are  recorded. 
Properly  Ccili'bratedj  the  Eppley  and  Kipp  pyranometers  have  an  accuracy 
eft  %, 

Photovoltaic  Instruments 

Ezperimenta  have  "been  performed  to  test  the  applicability  of  silicone 
photovoltaic  cells  for  solar  radiation  measurements  by  Tellott  (5)j 
Schoffer,   et,   al.   (6),  Banner^   et,  al.   (7)»   and  Whillier  (8).     A 
commercial  instrument,   the  Sol^A-Meter,   is  currently  manufactured  by 
Matrix,  Inoej  of  Phoenix,  under  a  license  from  John  Yellott. 

The  silicon  cell  consists  of  a  p-ci  junction  of  ultra-pure  silicon. 
Photons  striking  the  cell  cause  electrons  to  be  displaced,  which  is  a 
current.     Under  proper  loading  conditions,   the  current  flow  is  directly 
proportional  to  the  incident  radiation.     This  behavior  makes  the  silicon 
cell  a  candidate  for  solar  measurement. 

The  advantages  of  the  silicon  transducer  are  its  relatively  high 
electrical  output  and  its  long-tenn  stability  and  reliability.     Its 
disadvantages  include  a  slight  temperature  dependence  and  a  spectral 
sensitivity  which  departs  from  black  body  behavior.     Eitensive  experiments 
are  reported  in  the  references  above  which  compared  standard  Eppley  and 
Kipp  pyranometers  with  silicon  cell  pyranometers  (5j7}8),   These  experiments 
show  that  the  silicon  instruments  can  produce  results  within  t     5^.     Most 
of  the  inaccuracy  occurs  when  the  sun  is  at  low  angles.     The  variation  of 
calibration  with  optical  path  length  is  related  to  the  spectral  response 
of  the  cell.     This  effect  may  be  seen  in  the  calibration  data  for  the  SIMM" 
instruments  used  in  this  program. 
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Inrtrumenia  for  the  SI>M  P  ropy  am 

IrLB-trumenis  based  on  the  silicon  cell  transducer  were  chosen  for  iise  in 
the  SBM  prograia  for  two  reasonej     (l)   the  overall  cost  of  the  instruments 
Has  minimized  and  (2)   the  characteristicE  of  the  transducer  have  "been 
well  documented.     That  is,   the  silicon  cell  measurements  are  related  to 
the  standard  radiation  instruments. 

This  second  factor  is  very  important.     There  are  a  large  number  of 
photovoltaic  aad  photoresistive  transducers  available  at  the  present  time. 
Any  of  these  devicesf  when  connected  to  the  proper  circxiit  or  meter,  will 
give  some  change  in  output  with  variations  in  "light"  intensity.     The 
common  Cd-S  photographic  light  meter  is  an  example.     Due  to  linearity 
problems  and  spectral  response^   this  eystem  would  not  be  useful  for  enerer 
ir.eaGurements,  ^ 

The  immediate  purpose  of  the  SIMI.I  program  is  to  make  measurements  which 
will  accelerate  the  utilization  of  solar  energy^     It  is  therefore  important 
that  there  be  a  direct  relationship  of  the  measurements  made  to  the 
measurements  needed  for  solar  design^     Since  this  relationship  has  been 
documented  in  the  literature  ( 5^6, 7,8) ,  no  addJ.tional  development  or 
research  programs  are  needed  to  establish  the  "meaning*'  of  the  measurements. 

Prof.  John  Yellott  has  U.S.   patents  covering  some  aspects  of  solar 
measurin.g  instruments  using  silicon  cells.     I  have  described  the  SBM 
program  and  the  instruments  to  him.     He  has  assured  me  that  the  SBffl 
instruments  (.1^  :would  not  infringe  on  his  patents  and  (2)  would  be 
quixe  accurate  and  easy  to  calibrate,     . 

Description  of  the  SBM  Instrument  Fackap^e 

The  instrument  package  distributed  to  the  high  school  measuring  stations 
consists  of  tvjo  solar  radiation  measuring  units,  the  SBM-1   and  the  SBffl-IR. 
Photographs  of  the  instruments  are  shown  in  Figare     12. 

SIMM-I   is  a  hand-held  instrument.     It  consists  of  a  silicon  cell  mounted  on 
a  tiltable  plate.     The  plate  is  hinged  to  a  small  instrumentation  bo:c  which 
contains  the  indicating  meter.     The  instrument  is  self-powered  and  needs  no 
batteries.     It  is  used  in  this  program  as  a  back-^p  instrument  and  as  an 
educational  tool. 

The  SBOf-IR  instrument  incorporates  a  strip  chart  recorder  so  that  a 
continuous  record  of  solar  energy  can  be  obtained.     The  recorder  drive  is 
powered  by  100  vac«     The  transducer  is  separate  from  the  recorder  and 
consists  of  a  silicon  cell  mounted  on  an  angled  plate.     This  assembly  is 
designed  for  attachment  to  an  exterior  wall  or  roof.     The  recorder  is 
located  inside  the  building  on  a  convenient  shelf  or  cabinet. 

Transducer  Orientation  of  SDM~1R 

The  silicon  cell  in  the  SBM^IR  transducer  is  mounted  facing  south  and  at 
an  angle  of  60  degrees  from  the  horizontal.     This  orientation  is  chosen 
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because  it  is  nearly  ih.e  optinmm  winter  collection  angle  for  Montana* 
The  data  taken  at  this  angle  can  accxirately  lie  trasnfonned  to  other 
Birailar  angles.     Weather  liureaa  data  records  radiation  on  a  horizontal 
Burface.     Transforming  this  data  to  practical  collector  angles  produces 
error,     (in  Montana  in  the  winter  the  altitude  of  the  sun  is  about  20 
degrees.)     It  is  the  philosophy  of  this  program  to  cake  measurements 
which  will  implement  the  utilization  of  colar  energy.     This  goal  is 
achieved  by  taking  the  data  at  orientation  matching  the  ultimate  device. 


Errors  in  Data  Reduction 

The  new  direct  computer  data  reduction  has  improved  the  speed  and 
accxiracy  of  this  process.     VJe  have  tested  the  reproducability  of  the 
data  reduction  procedure  by  conducting  a  test  in  which  we  took  a  sample 
consisting  of  300  datardays  from  Billings  and  Bozeman.     This  300  day 
sample  was  read  two  times.     The  average  monthly  insolation  was  calculated 
sepatrately  for  each  run  and  then  compared.     The  variation  was  a  minimum 
of  less  than  1^  and  a  maximiim  of  ^.     The  average  variation  was  less  than 
2.5^. 
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Solar  Energy  Recording  System^   SIffl~1R 
20  units  currently  operating  across  Montana. 
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i,.       Manual  Solar  Energy  System,  SII-iil-l 
4   50  units  currently  located  across  Montana. 
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^PTIML  IMTER -STATION  SPACING 

The  concept  of  forecasting  time-series  events  is  not  new.   Almost 
any  variable  which  Is  observable  at  discrete  time  points  may  be  molded 
into  a  mathematical  model  which  may  closely  resemble  local-time-segment 
observations.   Once  a  valid  model  for  a  time  series  process  has  been  es- 
tablished, forecasting  analysis  may  render  useful  service  in  predicting 
performance  or  examining  system  modifications. 

Of  particular  interest  are  stationary  processes.   That  is,  those  pro- 
cessess  which  have  values  fluctuating  about  a  constant  mean.   These  processes 
are  unaffected  by  a  shift  of  the  time  origin,  and  may  be  modeled  by  taking 
observations  at  equal  time  periods.   Within  a  single  given  process,  a 
correlation  of  tliese  time  observations  (autocorrelation  may  appropriately 
be  used)  may  allow  construction  of  a  predictive  model  for  observations  at 
specific  time  increments  or  time  lags.   As  stated  by  Montgomery  and  Johnson: 
For  a  stationary  time  series  the  mean  is  just 

CO  oo 

t         J"0  3   t-j       j=0  J   t-3 

CO  CO 

and  since  the  sum  E.  ^=  A.  converges,  we  take  expectation  of  5!) .  _  ^^.e 

j=0  ^j       ■    '  ^  3=0  ^3  t-3 

CO 

term  by  term,  yielding  Z.  ^  ili.E(c   .)=0.   Thus  the  mean  of  the  station- 

3=0  \i    t-3 

ary  series  is         r./  \ 

E(x  )=m 

The  vj^riance  of  the  time  series  process  is   ^ 

Yo=V(x^)=E[x^-  E(x^)] 

«E[.Li0.e^  J^ 
J=0  3  t-3 

=E[.X„tJ».c   .   +  cross  products] 

3=0^3  t-3 

and  since  the  {e.}  are  independent,  the  cross  products  have  zero  expectation, 

X 

yielding 


oo      2 

The  variance  exists  only  if  5^.  r,=ilJ .  converges. 

3=0  J 

The  covariance  between  x  and  another  observation  separated  by  k  units 

t 

of  time  X  ,  is  called  "autocovariance"  at  lag  k,   nd  is  defined  as 
Yk-Cov(x^,  x^^^)=E[x^-E(x^)][x^^j^-E(x^^^)] 

*   From  Solar  Radia-tion  Data  Sources,  Applica-tionE  and  Network  Besign„, 
ERM  Contract  77-5-05-5362.  Carter,  et.al.,  197S. 


Tl.i.s  the  aL.-:^varlancG  is  just  like  the  covSance  of  two  random 
variables;  the  prefix  "auto"  merely  implies  that  we  are  referring 
to  the  covariance  of  any  two  observations  in  a  time  series  that  are 
k  time  periods  apart.   It  is  not  difficult  to  show  that  the  auto- 
covariance  at  lag  k  is 

2 
'k   e  j=o  J  j+k 

Within  the  framework  of  the  Box- Jenkins  methodology,  time  series 
models  are  characterized  by  their  autocorrelation  functions.   The 
correlation  between  two  random  variables,  say  W  and  Z,    is  defined 
as 

Cov(W,Z) 

p  ^^ . 

^'^Tvcw)  v(z)  ■ 

Thtis  the  autocorrelation  at  lag  k  refers  to  the  correlation  between 
any  two. observations  in  a  time  series  that  are  k  periods  apart.   That 
Is..   , 

Cov(x^  X  .,  )       Y^ 

t,  t+k'      _  »0 


1-2 


\ 


Vv(x  )  « 


^^\4k>     -Tl 


is  the  autocorrelation  at  lag  k.   A  graphical  display  of  v^   versus  lag 
k  is  called  the  autocorrelation  function  {y}  of  the  process..  Notice   ' 
that,  the  autocorrelation  function  is  dimensionless  and  that  -l<\i   <1. 
Futhennore,  Pj^=  P_j^ .  that  is,  the  autocorrelation  function  is  symmetric, 
so  that  it  is  necessary  to  consider  only  positive  lags.   In  general, 
when  observations  k  lags  apart  are  close  together  in  value,  we  would 
find  v^   close  to  1.0..  When  a  large  observation  at  time  t  i.s  followed 
by  a  small  observation  at  time  t+k,  we  would  find  jj  close  to  -1.0. 
If  there  is  little  relationship  between  observations  k  lags  apart,  we 
would  find  \}^   approximately  zero".'^ 

The  application  of  forecasting  technique  to  seasonal  actinometric 
data  by  L.  S.  Gandin  is  of  particular  interest  in  that  Gandin  computes  auto- 
correlation coefficients  of  observational  data  for  zero  time  lags.   These 
correlation  coefficients  are  then  graphically . displayed  as  a  function  of 


2 

mter-station  distance.    From  his  formula  for  interpolation  error  for 

square-grid-station  networks ,  he  then  computes  interpolation  error :^ 

^y^  (p/V2) 


E  =  1- 


1  +  2w  +  \}p^+n 


where: 


^n  =  the  measure  of  observational  error  (the  ratio 

of  the  mean  square  error  to  the  variance  of  in- 
terpolated magnitudes  here  taken  to  be  .25). 

p  =  the  correlation  coefficient  of  the  station  under 
scrutiny 

p  =  the  station-to-station  distance  or  the  area  attrib- 
uted to  a  single  station  within  a  square-grid  net- 
work. 

After  errors  have  been  computed  for  each  station,  they  are  plotted  as 
a  function, of  distance.   Optimal  inter-station  distance  may  then  be  deter- 
mined by  calculating  opti-aum  interpolation  error  for  the  number  of  states 
within  the  data  of  interest, 

*     1 

^     =  =  .25  for  n  =  2, 

2 
n 

and  reading  the  optimal  inter-station  distance  from  the  graph. ^ 

This  analysis  was  duplicated  for  four  actinometric  stations  in  South 
Carolina.   For  1970  seasonal  data,  individual  months  were  examined  ie.  Jan- 
uary, April  and  July.   Analysis  has  yielded  curves  which  are^^displayed, 
■^"-"  •   -   ■-   For  an  optimal  interpolation  error  of  .25,  the  inter-station 
distance  can  be  determined  from  Graph  2  to  be  approximately  75  kilometers. 

Examination  of  Graphs  )  and  Z   display  Gandin's  work  in  analyzing  A8 
stations  within  the  Soviet  Union.   For  an  optimal  interpolation  error  of 
.111,  Gandin  determines  the  inter-station  distance  to  be  150  kilometers. 
In  comparison  of  a  4-statiori  network  with  a  A8-station  network,  it 
must  be  said  that  only  a  demonstration  of  Gandin's  technique  has  been 
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Graph  S.-  Autocorrelation  function  of 

monthly  totals  of  global  radiation 

Summer. 
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Graph  ^  -   The  relationship  of  measure  of  optimum  Interpolation  error  in  the 
center  of  square.  Eg  and  of  triangle;,  G/i,  to  p  for  different  values 
of  the  measure  of  observation  error  r\.     Stations  are  supposed  to  be 
situated  at  the  apices  of  the  indicated  squares  and  triangles. 


Gandin's  curves  exactly  since  the  aerological  contours  of  the  regions  are 
different  and  since  the  period  of  records  are  not  the  sa.e.   Also  an  optical 
interpolation  error  value  of  .111  should  be  used  to  coincide  with  Gandin's 
network,  if  a  true  comparison  were  to  be  effective. 

In  conclusion,  to  further  investigate  this  technique,  data  fron,  a  larger 
nu..ber  of  stations  should  be  examined.   This  would  create  .ore  data  points 
for  curve-fitting  and  would  result  in  more  meaningful  results. 

The  purpose  of  explaining  this  procedure  is  to  provide  a  first  order 
approach  for  evaluating^ the  value  of  solar  radiation  measurements  interpolated 
to  intermediate  points.   This  will,  in  turn,  help  to  evaluate  or  Justify  the 
need  for  additional  data  in  the  network  area.  Additional  studies  of  this, 
and  other  methods  of  detennlning  optimum  spacing  of  stations  are  in  progress, 
and  the  reader  is  advised  to  seek  the  late.t  methods  for  determining  the 
quality  of  interpolated  data  between  stations. 
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roiogicai  oTi^i^^il^H^^TT— J7'    ^^^        ''/^°'=-    ^oTAFFendix  B.   World  hii^^^T" 
FTil.  ^^"    ^^^'   W^'^idMS^^^^^i^ii^Jr-Urganization   (1970), 
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Gandin,   loc  cit..    p.    2. 
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Introduction 

The  work  outlined  in  the  original  proposal  has  "been  accomplished.  The 
scope  of  the  work  was  increased  to  include  (a)  reduction  of  solar  data 
collected  betvjeen  grant  periods  and  (b)  preparation  and  reproduction  of 
a  summary  data  document.  This  was  done  with  no  increase  in  the  alloted 
"budget  amount. 

A  brief  ovejrview  of  the  scope  of  the  national  solar  radiation  data  base 
is  given  in  Figures  1  and  2.  Figure  1  shows  stations  for  which  reliable, 
global  radiation  is  archived  at  the  National  Climatic  Center.  Figure  2 
shows  current  stations  not  in  the  NOAA  network;  this  map  indicates  20 
stations  from  the  SIMM  network. 

Figure  2  shows  that  Montana,  Oregon  and  California  are  clearly  in  the 
lead  in  terms  of  solar  network  density,  Montana's  leadership  position 
is  due  to  the  support  of  SBIM  by  the  Alternative  Renewable  Energy  Sources 
Program. 

The  following  section  lists  milestones  accomplished  during  the  grant  year 
for  RAE  084^!LKS-772.  These  exhibits  are  intended  to  docxiraent  the  work 
done  on  this  grant. 

Milestones 

Exhibit  I;  Montana  Solar  Data  Manual 

The  Montana  Solar  Data  Manual  was  compiled  in  late  spring  and  early  summer 
during  an  extension  on  the  grant.  This  170  page  Manual  summarized  solar 
data  collected  during  the  past  two  grant  years  and  includes  auxiliaiy 
information  useful  to  the  solar  designer. 

A  total  of  350  copies  of  this  Manual  were  reproduced:  250  copies  were 
sent  to  BNR,  50  copies  were  distributed  to  participating  high  schools 
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SOLAR  RADIATION  STATIONS  WITH  DATA  ARCHIVED  AT  THE  NATIONAL  CLIMATIC 
CENTER,  ASHEVILLE,  NORTH  CAROLINA.   NUMBERS  INDICATE  MORE  THAN  ONE 
STATION  IN  AREA. 
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FIGURE  2  -  SOLAR  RADJATION  STATIONS  NOT  IN  THE  BASIC  NATIONAL  WEATHER  SERVICE, 
NOAA,  NETWORK.   ASTERISKS  DENOTE  DOE  RESEARCH  AND  TRAINING  SITES.   NUMBERS  INDICATE 

MORE  THAN  ONE  STATION  IN  THE  AREA. 


and  40  copies  were  distributed  primarily  to  professionals  and  government 
agencies  outside  of  Montana. 

(Notes  Copies  of  the  Manual  are  not  attached  to  this  report  as  tbey  are 
available  at  IWR.) 

Exhibit  II:     The  Montana  Solar  Radiation  Measurement  Network     (attached) 

This  paper  was  presented  at  Solar   »79  Northwest   in  Seattle  on  August   11 . 
This  meeting  was   jointly  sponsored  "by  the  U.S.   Department  of  Energy, 
Bonneville  Power  Adjninistration,  Vfashington  State  Energy  Office  and 
Oregon  Department  of  Energy.     Over  100  technical  papers  were  presented 
for  the  (estimated)    1000  people  attending  the  conference. 

The  paper  describes  the  SBM  program  operation,    instrumentation  and 
data  output. 

E3chibit  III;     Handbook  of  Solar  Energy  Experiments  (attached) 

This  Handbook     was  prepared  as  an  integral  part  of  the  SIMM  program  for 
distribution  to  the  participating  high  schools.      The  earlier  version  of 
this  Handbook  has  been  considerably  enlarged  this  year  by  including 
several  new  solar  energy  experiments.     The  new  experiments  focus  on 
solar  utilization.     They  include  instructions  for  building  and  testing 
simple  active,   passive  and  concentrating  collectors.     Most   experiments 
may  be  carried  out  as  simple  demonstrations  or  as  more  detailed  tests. 
The  solar  basics  are  emphasized  to   encourage  the  students  to  grasp  the 
concepts  of  solar  utilization.     The  30  schools  having  insolation  recorders, 
as  well  as  the  schools  having  manual  instruments,   receive  the  Handbook. 

Exhibit  IV J      "Solar  Radiation  Observation  Stations,  Updated  to   1979"   (attached) 

This  document  was  prepared  by  the  Johnson  Environmental  and  Energy  Center 
of  the  University  of  Alabama.     The  purpose  of  this  effort  was  to  catalog 
stations  within  the  U.   S.     where  solar  insolation  was  being  measured.     The 
maps  in  Figures  1    and  2  were  taken  from  this  document  and  all  thirty  of 
the  stations   in  the  SIMM  network  are  listed. 


Exhibit  V;     NOAA  Calibration  (attached)  ' 

The  silicon  cell  transducer  used  in  the  SBIM  program  was  calibrated 
by  Edwin  Flowers  of  NOAA   for  approximately  two  weeks  against  the  NOAA 
reference  standard.     This  data  will  allow  us  to  more  accurately  establish 
the  accuracy  of  our  data  in  the  future.     As  a  parallel   effort,    the  Kipp 
pyranometer  which  had  been  used  as  our  secondary  standard  was  recently 
re-calibrated  by  NOAA.     The  calibration  of  the  Kipp  had  not  changed 
significantly.     These  efforts  in  conjunction  with  our  clear-day  procedure 
insure  the  integrity  of  our  data. 

Conclusion  . 

In  addition  to  these  major  milestones  there  has  been  a  gradual  increase 
in  awareness  of  the  SIMM  program  in  people  outside  Montana.   I  have  made 
two  presentations  in  VJashington  arid  Oregon  to  regional  groups  involved 
in  radiation  measurement.  We  have  also  answered  numerous  requests  for 
information  from  individuals  within  DOE,  SEHI,  State  Energy  Offices  and 
Universities  concerning  the  SUM   program.  I 

Measurements  and  data  compilation  are  currently  being  continued  at  the 
30  stations.  During  the  summer  a  technician  visited  all  stations,     : 
checked  the  instruments  and  made  necessary  repairs.  During  the  coming 
year  more  work  will  be  done  on  calibration  and  ideas  for  improving  the 
Data  Manual  will  be  solicited. 

Budget 

Table  I  outlines  the  expenditures  as  of  March  31?  1979»  At  this  point 
there  was  a  balance  of  $3,952.48  and  the  goals  of  the  proposal  had  been 
accomplished  with  the  exception  of- the  data  handbook.  At  this  time  a 
request  v/as  made  to  expand  the  scope  of  the  proposal,  using  these  funds 
to  pay  for  the  extra  work.  This  recfuest  was  approved  by  DNRC. 

There  were  t^^:o  main  tasks  in  the  expansion.  The  first  task  v;as  to 
reduce  strip  chart  data  which  Fowlkes  Engineering  had  collected  during 
late  1977.  This  data  had  never  been  reduced  because  it  was  accvimulated 


during  a  gap  in  the  funding  between  the  first  and  second  year  of  the  grant. 

The  second  goal  was  to   expand  the  proposed  50  copies  of  the   I978  data 
document  into  350  copies  of  a  much  broader  document,    the  Montana  Solar 
Data  Manual,    containing  both  1977  a^^  1978  data  along  with  auxiliary 
design  information. 

Table  II   itemizes  the  use  of  the  $3952.48  for  accomplishing  these 
expansions.     It  vjill  be  noted  that  the  expansion  cost  exceeded  the  money 
available  by  $663,77,    ^'^^  primarily  to  a  large  amount  of  computer  time  to 
output  the  Data  Manual,     This  overrun  will  be  absorbed  by  Fowlkes 
Engineering. 

Final  Budget  Request 

Enclosed  is  a  grant  request  for  $3547.50.     This  is  the  amount  remaining 
in  the  grant  after  the  previous  quarterly  advances.      The  final  budget 
distribution  was  much  the  same  as  the  original.     The  largest  transfer 
was  $673.98  from  maintenance  into  Data  Manual   reproduction:   the  remaining 
transfers  were  typically  around  $200. 


THE  MONTANA  SOLAR  RADIATION  HEASURIHG  NETHQItK 
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ABSTRACT 

This  paper  doocrlbaB  the  operation  and  results  of  a  current 
state-supported  program,   Solar  Insolation  Measurement,  Mon- 
tana (SOT!).     The  goal  of  this  program  is  to  measure  and 
document  the  solar  insolation  available  at   JO  locations 
within  the  state  of  Montana.     This  solar  insolation  data 
base  will  be  useful  to  designers  of  solar  enorgjr  conversion 
systems.     The  existence  of  this  data  base  is  intended  to 
stimulate  the  utilisation  of  solar  energy  in  Montana.     To 
date,   over  20,000  solar  insolation  data-days  have  been  col- 
lected.    A  sujwnaiy  of  most   of  this  data  is  available  to   the 
public  in  the  recently  published  Montana  Solar  Data  Manual, 
rig.    1    (Fowllces,    1979).  """" ~ 
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Fig.  1.  Cover. of  Montana  Solar  Data  Manualf  illus- 
trating solar  insoiation  records  (including  eclipse 
on  Feb.    26,    I979). 


INTRODUCTION 

Early   in   1977  the  SIMM   program  (Solar  Insolation  Measure- 
ment, Montana)   was  initiated  by  Fowlkes  Engineering  under  a 
p-ant   from  the  Alternative  Renewable  Energy  Sources  Program 
of  the  Montana  Department  of  Natural  Resources  and  Conser- 
vation.     Twenty   instruments  to  measure  solar  radiation  were 
sited  across  Montana  during  1977  and  an  additional  ten 
instruments  were  sited  early  in   I978.      In  addition  to  the 
thirty  stations  having  continuous  insolation  recorders 
there  are  currently  twenty-five  stations  -which  report  the 
readings  of  manual    instruments,  Fig.    2. 
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•     Continuous  Data  Recording  Station 
o     Manual  Data  Station' 
Pig. 2.     Montana  Solar  Insolation  Measurement  Network 

All  of  the  instruments  are  located  in  public  high  schools. 
At   each  school  there   is  a  volunteer,  usually  a.  science 
teacher,  who  is  responsible  for  the  instrument.     The  vol- 
unteer periodically  replaces  the  chart  paper  in  the  record- 
er and  mails  the  completed  strip  charts  to  our  office  for 
processing.     The  total  time  requll-ad  for  the  volunteer  to 
do  these  tasks  usually  amounts  to   less  than  five  hours 
each  year. 

The  solar  insolation   strip  chart   reoordors  are  often  loca- 
ted in  a  laboratory,   office,   or  olasaroom  where  they  can 
be  seen  by  the  students.     Some  teachers  make  specific  use 
of  the  insolation   records  in  their  olassss. 

All  high  aohools  participating  in  ths  SIMM  program  pooaive 
a  copy  of  The  Handbook  of  Solar  F.nmer  E:toerJmfmtfl  and  a 
manual   solar  insolation  meter.     This  Handbook  was  written 
by  the  author  of  this  paper  as  an   integral  part  of  the  in- 
solation measurement  program.  The  80  page  Handbook  contains 
a  review  of  solar  basics  and  ten  outlined  experiments 
ranging  from  insolation  measurement  to  testing  a  solar 
collector. 
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INSOLATION  ON  A  TILTKU  SURFACE 

Meteorological  insolation  data  Id  uoually  recorded  on  a 
horizontal  surface  due,  in  part,  to  the  fact  that  the  solar 
collector  of  interest  to  meteorolofjiots  ia  the  horizontal 
surface  of  the  earth.  For  northern  U.S.  latitudes  the 
solar  designer  vdll  normally  tilt  the  collector  plane  be- 
tween 45*^  and  90*^  up  from  the  horizontal  to  increase  the 
average  radiant  flux  density  intercepted  by  the  collector. 
Procedures  exist  for  estimatin^^  the  inoolation  on  a  tilted 
surface  given  the  measured  insolation  on  a  horizontal  cur- 
face.  However,  the  net  uncertainties  of  these  estimates 
in  northern  H.S.  latitudes  during  the  winter  will  be  around 
■^  20^.  Some  uncertainties  arise  from  the  pyranonieter 
errors  at  low  solar  altitudes.   Additional  uncertainties 
exist  due  to  reflected  radiation  from  snow  covered  fjround 
which  will  be  intercepted  by  steep  solar  collectors  but  will 
not  be  intercepted  by  a  horizontal  radiation  transducer. 

Since  the  goal  of  the  SIMM  pro-am  is  to  develop  a  solar 
insolation  data  base  for  designin,^  solar  conversion  devices 
we  measure  insolation  on  a  surface  tilted  60°  up  from  hori- 
zontal and  facing;  due  couth.  The  transducer  hence  becomes 
a  model  of  a  typical  solar  collector  which  might  be  used  at 
thi-«  latitude  to  collect  heat  in  the  winter-,  The  resultin/j 
60°  solar  radiation  data  can  often  be  used  directly  or  can 
be  transformed  to  neighboring  tilt  angles  with  relatively 
small  increases  in  uncertainty.  The  effects  of  snow  and 
dirt  on   a  collector  are  approximately  modeled  by  this  trans- 
ducer orientation  and  arc  an  integral  part  of  the  insolation 
data. 
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low  energy  density  on  cloudy  days.  Similarly,  amplitude 
errors  near  eunrlae  and  sunset  have  a  small  effect  on  daily 
total  energy  density:  less  than  2^  on  a  clear  day. 

The  published  absolute  accuracy  of  the  atrip  chart  recorders 
is  y^.     This  error  is  reduced  significantly  since  the  re- 
corders are  directly  calibrated  in  conjunction  with  the 
transducers.  The  readability  of  the  recorder  during;  cali- 
bration, using  a  magnifier,  is  better  than  -  l^j. 

The  electrical  circuit  used  is  sensitive  to  changes  in  lead 
resistance.  This  factor  has  probably  been  the  primary 
source  of  error  in  this  measurement  program.  Resistance 
variations  are  usually  due  to  improperly  soldered  lead  wire 
connections.  The  {jpi^adually  increasing  contact  resistance 
leads  to  a  gradual  reduction  in  sensitivity  of  the  instru- 
ment.  In  a  few  cases  this  problem  has  been  responsible  for 
temporary  errors  of  10  and  205*.  Vie  are  currently  inspect- 
ing and  servicing  the  instruments  at  all  30  stations. 


DATA  REDUCTION 

Strip  chart  data  is  processed  by  tracing  the  curves  on  a 
digiti7,er  platen  connected  to  a  microcomputer,  Fig.  4.  The 
resolution  of  the  platen  is  0.1  mm  and  each  curve  is  repre- 
sentad  by  200  to  400  coordinates.   The  reading  error  due  to 
resolution  is  thus  less  than  1^,  The  primary  source  of 
data  reduction  error  is  the  repeatability  of  the  curve 
tracing  itself.  On  clear  daySf  where  the  curve  is  a  well 
defined  line,  the  reproducibility  averages  1^.  On  cloudy 
days  the  reproducibility  can  reach  5  *o  20^.  This  is  due 
to  the  fast  response  of  the  silicon  cell  which  producee  a 
very  rough,  jagged  curve  on  some  cloudy  days,  Fi^:.  1, 


The  solar  insolation  is  sensed  by  a  silicon  cell  mounted  in 
an  8  cm  square  aluminum  plate  and  covered  with  glanSfl  shown 
in  Pig.  3.  The  transducer  is  attached  to  a  bracket  outside 
the  building  and  a  wire  runs  to  a  strip  chart  recorder  loca- 
ted inside  the  building,  V/e  fabricated  all  the  transducers 
used  in  this  program. 


Fig.  4. 


Computer  system  used  to   read  strip-charts  and 
process  solar  insolation  data 


Fig. 3-    Instrumentation   for  recording  solar   insolation 


Silicon   coll   characteristico   relating  to    solar   fliDc  measure- 
ment are  well  documented  in  the   literature  and  several  com- 
mercial  silicon  cell   instruments  are  available  (Selcuk, I962). 
Briefly,   the  silicon  cell   characteristics  causing  departures 
from  pyranometer  readings  are  due  to   spectral   response, 
temperature  changes  and  cosine  response.   The  epeotral   re- 
sponse will  cause  the  silicon  cell  to   road  measurably  hi^or 
on   cloudy  days  and  within   2  hours  of  sunrise  or  sm;set.   The 
instruments  used   in   this  program,    for  example,   will   read 
about  0.015  kiv/m     high  under  conditions  of  80^  beam  attenu- 
ation due  to   clouds.   The  amplitude   error  in  this  condition 
is  ±.  6^,  The  effect  of  this  spectral   error  on  the  monthly 
average,   daily  total   energy  density  is  typically  reduced  to 
less  than   3/^.     This   reduced  error  is  due  to  the   relatively 


During  the   last  year  reproducibility  tests  have  been   run 
on   over  8OO  data  days   (4?^  of  our   total   data  base).   The   re- 
producibility of  the  monthly  average,    daily  total    energ/ 
density  averages  better  than  -  2fo,     These  monthly  values 
often    reproduce  within    1^  due   to  the   compensation  of  posi- 
tive and  negative   errors  when   averaged. 
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The  data  are  stored  on  hard-eeotored,    5^  inch  maEnetio  disk- 
ettes using  the  North  Star  Disc  Operating  System  (a  popular 
microcomputer  format).     Each  single-density  diskette  con- 
tains 3   station-years  of  solar  data.      This  data  has  been 
used  by  Fowlkes  Engineering  as  input  to   solar  heating  simu- 
lation models.      Copies  of  these  diskettes  can  be  supplied 
at  nominal  cost   for  users  wanting  to   implement  their  own 
simulation  programs. 


The  clear  day   procedure   is  used  to   evaluate  all  of  our  data. 
If  the   experimental   clear  day  values  are  consistently  in 
disagreement  the  data  is  scaled  by  on  appropriate  correction 
factor.      An   exajnple  of  this  is  given  in  Pig.    5b  for 
Glendive,  Montana.      In  this  case  the  experimental   clear  day 
values  appeared  uniformly  low  by  4-5^  and  all  the  Olendive 
data  for  that  year  was  therefore  scaled  up  by  4. 3^.     The 
net  effect  of  this  procedure  is  to  insure  that  large  errors 
do  not  occur  duo  to   instrumentation  calibration  drift. 


DATA  ACCimACY 

All  transducers  are  calibrated  in  conjunction  with  the  strip 
chart  recorder  by  comparing  the  SBIM    instrument  to   a  Kipp  or 
Eppley  pyranometer  on  a  clear  day  near  noon.     The  calibra- 
tion factor  for  each  instrument   is  then  determined  from  this 
comparison.     Although  this  5)rocedure  established  an  initial 
calibration  reference  point   it  does  not  account   for  the 
effects  of  drift,    aging,    etc.   which  may  occur  on  the  site 
after  the   instrument   is  installed.      For  long-term  control  of 
data  accuracy  we  make  use  of  a  clear  day  technicfue.     This 
technique  parallels  the  methods  used  by  D.O.E.    in  their 
recent   rehabilitation   effort  to  create  the  SOLMET  data  base 
(DOE,    1979). 

The  essence  of  the  clear  day  procedure  is  to  compare  the 
measured  clear  day  radiation  to  the   radiation  predicted  by 
a  theoretical  clear  day  model   and  then  make  appropriate 
corrections.     VJe   studied  two  clear  day  models.     The  first 
model  was  provided  by  Doug  Hoyt  of  NOAA   (Hoyt,1978).  This 
model  was  used  in  the  SOLMET  data  rehabilitation.     The 
second  model  was  based  on  the  ASHRAE  procedure  and  constants 
with  a  clearness  modification  to  account   for  altitude. 

As  part  of  our  data  reduction  procedure  we  pick  off  the  max— 
irarnn  radiation  at  solar  noon  on  all  clea.r  days  and  store 
this  value  on  the  computer.      The  ASHRAE  clear  day  calcula- 
tions were  also  programmed  on  our  computer.     Doug  Hoyt  pro- 
vided us  with  clear  day  noon  maximums   from  his  model   for  a 
60°  tilted  surface.     A  graph  of  these  two  clear  day  models 
along  with  experimental   cleax  day  values  are   shown   in  Vig.  5a 
for  Great  Falls,  Montana.     Based  on  this   comparison  we  con- 
cluded that  the  modified  ASHRAE  clear  day  model  could  be 
used  to  maintain  the  integrity  of  our  data. 
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(b)  Glendive  data  for  1978 
illustrating  use  of  scale  factor 


"320       360 


Fig.    5.     Comparison  of  clear  day  experimental  and 
tneoretical  values 


DATA  OUTPUT 

The  solar  insolation  data  output  from  this  program  has  been 
summarized  in  the  Montana  Solar  Data  Manual  which  will  be 
distributed  to  libraries,  extension  sorvices,  engineers, 
architects  and  solar  designers  this  summer.  This  Manual 
includes  the  two  years  of  solar  insolation  data  collected 
to  data  for  all  the  stations  in  the  network.  Data  on  the 
frequency  of  occurence  of  cloudy  days  has  also  been  inclu- 
ded to  aid  in  the  sizing  of  heat  storage  Bystems.  Average 
climatological  data  from  National  Weather  Service  records 
is  reproduced  in  the  Manual  for  the  convenience  of  the 
designer.  Also  included  in  the  Manual  is  a  brief  review  of 
solar  design  methods,  economic  analysis,  a  list  of  refer- 
ences, conversion  tables  and  a  diecuBsion  of  the  accuracy 
of  the  solar  data. 

The  solar  data  is  presented  in  tables  in  the  form  of  weekly 
and  monthly  averages  of  the  daily  total  energy  density  for 
each  station.  Also  given  is  the  average  intensity  spectrum 
of  the  daily  total.  This  is  prosentad  as  a  percentage  of 
the  daily  total  energy  density  which  occured  at  radiation 
intensities  above  12  threshold  values  ranging  from  0.1  to 
1.2  kH/m  .  Due  to  losses,  solar  collectors  require  a 
radiation  of  some  minimum  intensity  before  they  can  produce 
positive  output.  This  threshold  may  be  from  0.1  to  0.4 

SOLA}!  KADIATIOH  DATA  fori  HELliMA,  MONTANA  1978 
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Fig.    6.     Sample   solar  data  from  manual 
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IcH/m^  depending  on  the  collector  deoign  and  ambient  tempera- 
ture.    Given  the  threshold  for  the  collector  the  correspond- 
ing percentage  may  be  read  from  the  spectrum  table.     This 
percentage  multiplied  by  the  total   energy  density  will  rep- 
resent the  maximum  usable  solar  energy  available  to   the 
collector. 

Figure  6  reproduces  a  page  of  weekly  average  data  from  the 
Vlonta-ia  Solar  Pata  Manual   for  Helena,  Montana  during  197H. 
The   first   column  lists  the   starting  date   for  the  week  and 
the  second  column  lists  the  average  daily  total  energy  den- 
sity,  kW-hr/m^-day,    during  that  week.     The  next  twelve  col- 
umns list   the  averaga   spectrum  during  the  week.     Similarly, 
more  compact  tables  listing  monthly  averages  are  also   inclu- 
ded in  the  Manual.      Parentheses  around  the  daily  averages 
denote  that  at   least  one  day  of  data  was  missing  durxng  the 
averaging  period. 

The  current  Manual   contains  solar  data  from  the   first  two 
years  of  the  measurement   program,    1977  and  1978.     These 
measurements  are  being  continued  throueJi  the  next  year  and 
the  Data  Manual  will  be  updated  at  that  time.     Assuming 
funds  are  available,    this  program  will  be  continued  to   a 
point  where   long-term  average  insolation  data  con  be 
reliably  documented. 


RELATIVE  SOLAR  POTENTIAL 

The  solar  data  for  the   last  two  years  has  been  summarized  to 
give  a  preliminary  overview  of  the  relative  potential   for 
solar  utilization   in  Montana.     Solar   insolation  data  and 
long-term  average  degree-day  data  were  averaged  for  the 
8  month  heating  season,   October  through  May.     The  average 
solar  insolation  was  multiplied  by   on    efficiency   factor  and 
then  divided  by  the  total  heating  season   degree-days.     The 
efficiency   factor  was  derived  for  each   station  based  on  a 
nominal  collector  efficiency  curve.      These   efficiency  values 
were   similar  for  all   stations.      In  areas  with  high  relative 
solar  potentials,    a  given  collector  design  would  supply  a 
hi^er  percentage  of  the  annual  heat   load.      Further  refine- 
ment of  the  solar  potential  of  different  areas  of  Montana 
will  require  both  more  detailed  models  and  a  largsr  solar 
insolation  data  base.     The  current   results  are  presented  in 
Fig.    7. 


ACKNaiLEDGIWEWTS 

This  measurement  program  is   funded  by  a  grant   from  the 
Alternative  Renewable  Energy  Sources  Program  of  the  Montana 
Department  of  Natural  ResouroeB  and  Conservation.     The 
support  and  advice  of  Dana  Ounderson,   Administrative  Officer, 
is  acknowledged.     The  assistance  of  the  volunteers  at  the 
measurement   stations   is  crucial  to   the  buocosb  of  this  pro- 
gram and  their  contributions  are   gratefully  acknowledged. 

Several  members  of  the   staff  at  Fowlkes  Engineering  have 
contributed  to  this  solar  measurement   program  and  this 
paper.     Pat  Sullivan  has  carefully  logged  and  organized 
strip  chart   records  for  over  15,000  data  days.      Pat  also 
typed  this  paper  and  made  editorial  corrections.     Sharon 
Kleingartnar  is  acknowledged  for  her  important   role  in 
entering  the  data  into  the  computer  for  processing.  Kellee 
Crisafulli   and  Ross  Nelson  have  written  numerous  compruter 
programs  to  process  and  verify  the   solar  data. 


RliFIJREHCES 

Fowlkes,   CharleseW.,  Montana  Solar  Data  Manual,    available 
from  Montana  Dept.    of  iTatural  Resources  i  Conservation, 
Alternative  Renewable  Energy  Sources  Program,    32  S.   Swing, 
Helena,  MT   5971 5>    1979- 

Hoyt,   Douglas  V.,    "A  Model   for  the  Calculation  of  Solar 
Global  Insolation",  Solar  Energy,  Vol.    21.  No.    1    (l978) 

Selouk,  K.    and  Tellott,   J.    I.,    "Measurement  of  Direct,  Dif- 
fuse,   and  Total  Solar  Radiation  with  Silicon  Photovoltaic 
Cells",  Solar  Energy,  V(?l,    6,  No.    4  (l$62) 

U.   S.   Dept.    of  Energy^iv.    of  Solar  Technology,    "Hourly 
Solar  Radiation-Surface  Meteorological  Observations", 
SOLMET     Final  Report,  Vol.    1   (1979) 


I..O.40-0.46,   11-0,35-0.39,   III-O.30-0.34,   IV.0.20-O.29 
Fig. 7.     Map  showing  four  regions  of  relative   solar/load  ratio 


266 


laniBIT  W:      SOLAR  RADIATION  OBSERVATION  STATION! 
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FAT ION  LOCATION 
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Will" 18' 
1015  M 
FE 

Hamilton 
N  i{6°15' 
Wlli{'09' 
1097  M 
FE 


Pyranometer  (T) 
Silicon  -  SIMM  -  IR 


Pyranometer  (T) 
Silicon  -  SIMM  -  IR 


Strip  chart 
Rustrack  Recorder 
(15) 


1977  - 


Strip  chart       1977 

Rustrack  Recorder 

(15) 


Harlowton 
N  i}6'26' 
W109'50' 
1270  M 
FE 


Pyranometer  (T) 
Silicon  -  SIMM  -  IR 


Strip  chart 
Rustrack  Recorder 
(15) 


1977 


Havre 
N  ^8-3V 
,,-4/109-40' 
38  M 
.£ 


Pyranometer  (T) 
Silicon  -  SIMM  -  IR 


Strip  chart 
Rustrack  Recorder 
(15) 


1977  - 


STATION  LOCATION 


Helena 
N  46-35' 
W112°02' 
1257  M 
FE 


RADIOf^ETER 

RECORDER 

PERIOD  OF 

TYPE,  MFG.  MODEL 

RECORD  FORM 

RECORD 

MONTANA 

(CONTD) 

Pyranometer  (T) 

Strip  chart 

1977  - 

Silicon  -  SIMM  -  IR 

RusTRACK  Recorder 
(15) 

Jordan 
H  47-19' 
W105'54' 
853  M 
FE 


Pyranometer  (T) 
Silicon  -  SIMM  -  IR 


Strip  chart 
Rustrack  Recorder 
(15) 


1977 


Kali  spell 
n  48-13' 
'//114"i5' 
902  M 
FE 


PYRANOMEtER  (T) 

Silicon  -  SIW  -  IR 


Strip  chart 

Rustrack  Recorder 
(15) 


1977 


LEW  IS TOWN 

N  47-04' 
W109-25' 
1207  .1 
re 


Pyranometer  (T) 
Silicon  -  Sim  -  IR 


Strip  chart 
Rustrack  Recorder 
(15) 


1977 


Ll3BY 

N  48-24' 
W115'35' 
625  M 


Ft 


Pyranometer  (T) 
Silicon  -  SIMM 


Strip  chart 
IR     Rustrack  Recorder 
(15) 


1977  - 


Livingston 
N  45-39' 
WHO- 34' 
1369  M 
FE 


Pyranometer  (T) 
Silicon  -•  SIMil  -  IR 


Strip  chart 
Rustrack  Recorder 
(15) 


1977 


Miles  City 
N  46-25' 
W105'49' 
72  M 
FE 


Pyranometer  (T) 
Silicon  -  SIMM 


Strip  chart 
IR     Rustrack  Recorder 
(15) 


1977  - 


TION  LOCATION 


RADIOr.ETER 
TYPE.  MFG.  MODEL 


RECORDER 
RECORD  FORM 


PERIOD  OF 
RECORD 


Missoula 

982  M 

FE 

Plentywood 
U  ^8-47' 
W10^°3ii' 
517  M 
FE 


mam 

(CONTD) 

Pyranometer  (T) 
Silicon  -  SIMJ^  -  IR 


Pyranometer  (T) 
Silicon  -  SIM.f^  -  IR 


Strip  chart       1977 

RusTRAC.K  Recorder 

(15) 


Strip  chart       1977 

Rustrack  Recorder 

(15) 


POLSON 

H  il7°41' 
WIWIO' 
899  M 


Pyranometer  (T) 
Silicon  -  SIM^  -  IR 


Strip  chart 
Rustrack  Recorder 
(15) 


1977  - 


i\£D  Lodge 
H  A5°ll' 
W109°15' 
1691  M 
FE 

Sidney 

ti  ^7'^2' 

588  M 

FE 


Pyranometer  (T) 
Silicon  -  SIMJ1  -  IR 


Pyranometer  (T) 
Silicon  -  Sim  •   IR 


Strip  chart       1977 

Rustrack  Recorder 

(15) 


Strip  chart       1977 

Rustrack  Recorder 

(15) 


Thompson  Falls 
N  «<7'35' 
m5-20' 
751  M 
FE 

West  Yellowstone 

Wlll-04' 

■31  M 


Pyranometer  (T) 
Silicon  -  Sim  -   IR 


Pyranometer  (T) 
Silicon  -  SIMM  -  IR 


Strip  chart       1977  - 

Rustrack  Recorder 

(15) 


Strip  chart       1977  - 

Rustrack  Recorder 

(15) 


U.S.  DEPARTMENT  OF  COMMERCE 

SUational  Oceanic  and  Atmospheric  Administration 

ENVIRONMENTAL  RESEARCH  LABORATORIES 
Boulder.  Colorado    B0302 


EXHIBIT  V:      NOAA  CALI15KATI0N 


RADIATION  INSTRUMENT  CALIBRATION 


SERIAL  NO. 


(None) 


INSTRUMENT  TYPE  Pyranometer 
MANUFACTURER     Fowlkes  Eng. 


MODEL  •  Photocell 

THIS  INSTRUMENT  HAS  BEEN  COMPARED  WITH  A  NOAA  REFERENCE  STANDARD 

THE  DERIVED  CONSTANT  IS: 
80,3    X  10"^  volts/watt  meter~2 
56^0      millivolts/langley  niin"^ 


ON  THE  ABSOLUTE  RADIATION  SCALE  1975 


April  10-19,  1979 


DATE  OF  TEST 


Edwin  C,  Flowers 
SOLAR  RADIATION  FACILITY 


Calibration  performed  for: 

Char less  W.  Fowlkes 
31  Gardner  Park  Drive 
Bozeman,  MI  59715 


%^ 


^ 
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U.S.  DEPARTIVlEr\IT  OF  COiVlMERCE 

rUational  Oceanic  and  Atmospheric  Administration 

ENVIRONMENTAL  RESL-ARCH  lABORATORIES 
Boulder,  Colorado  B03U2 

January  1,  1977 


Calibrations  performed  by  the  Solar  Radiation  Facility  are  on  the 
Absolute  Radiation  Scale„   This  scale  is  embodied  in  the  instrument 
PACRAD  III,  an  active  cavity  radiometer  maintained  by  the  World 
Radiation  Center,  Davos,  Switzerland.   (For  reference  see:  Pyrheliometer 
Comparisons  Davos,  1975,  Results,  Working  Report  No.  58,  Swiss  Meteoro- 
logical Institute  Zurich,  Switzerland,  1976.) 

The  relationship  between  calibration  factors  based  on  the  Absolute 
Radiation  Scale  (ABS)  and  factors  based  on  the  International  Pyrhelio- 
metric  Scale  (IPS)  as  maintained  by  NOAA,  is: 

C(ABS)  =  0-975  C(ips) 
This  means  that  measurements  made  with  instruments  calibrated  by  NOAA 
on  the  IPS  will  yield  irradiance  values  approximately  2.57o  lower  than 
measurements  made  with  instruments  calibrated  by  NOAA  on  the  Absolute 
Scale.  N 


Edwin  Flowers,  Chief 
Solar  Radiation  Facility 


-J^AY^ 


FINAL  REPORT 
TO 
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Effective  Date:  4/17/79 


from 

Fowlkes  Engineering 

31  Gardner  Park  Drive 

Bozeman,  MT       59715 


July  1980 
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1»0  IHTROmCTION 

The  emphasis  diiring  the  third  year  of  the  Solar  Insolation  Measure- 
ment, Montana  (SIMM)  prqgram  has  "been  to  (a)  continue  the  solar  data 
collection  processg  (b)  make  the  current  data  avedlable  to  the  public 
and  (c)  document  the  quality  and  accuracy  of  the  date*  Specific  mile- 
stones and  exhibits  are  included  in  this  report  to  demonstrate  these 
accomplishments. 


2.0  BUILDING  THE  DATA  BASE 

During  the  grant  year  over  9»000  charts  of  daily  solar  radiation 
were  collected  from  thirty  stations  in  the  field.  These  charts  were 
checked  for  errors,  analyzed^  read  onto  a  computer  system  and  filed  for 
future  reference.  The  computer  data  base  was  subsequently  checked  for 
errors  and  archive  copies  of  the  data  were  made  and  stored. 

This  data  was  output  in  the  fonn  of  tables  of  solar  radiation  and 
cloud  cover  at  each  of  the  thirty  sites.  The  tables  were  reproduced  and 
bound  to  form  the  1979  Addendum  to  the  Montana  Solar  Data  Manual.  Three 
hundred  copies  of  the  Addendum  were  delivered  to  the  Department  of 
Natural  Resources  and  Conservation  as  part  of  the  grant  agreement.  Ex- 
hibit I  is  a  copy  of  this  document. 


3*0  INSTRUMENT  MABPTHTAUCE 

Throu^out  the  year  the  instruments  in  the  field  were  repaired  and 
maintained.  Some  of  the  recorders  have  run  continuously  for  three  years 
at  this  point  and  occasioneilly  motors  will  bum  out  or  the  chart  drive 
linksige  will  fail.  Part  of  the  maintenance  program  includes  routijie 
phone  calls  to  the  volunteers.  During  the  summer  of  1979  a  technician 
visited  all  30  sites  in  person.  He  performed  routine  maintenance  and 
checked  the  orientation  and  calibration  of  the  solar  recorders.  Details 
of  the  installation  were  documented  on  the  form  shown  in  Exhibit  II. 


4.0  TECHNICAL  PRESETTTATIONS 

Teclmical  papers  on  this  solar  measurement  program  were  presented 
and  published  at  two  national  conferences:  Solar  Northwest  79  in  Seattle 
and  the  annual  International  Solar  Energy  Society  meeting  in  Phoenix  in 
1980,  These  presentations  have  allowed  a  wide  audience  of  solar  experts 
to  review  our  methods  and  results.  Their  response  has  heen  favorable, 
A  reprint  of  the  most  recent  paper  is  given  in  Exhibit  III.  This  paper 
provides  a  goodj  brief  description  of  the  program. 


5.0  SOLAR  DATA  mmJAL   DISTRIBUTIOM" 

During  the  past  year  the  Manual  has  been  sent  to  individuals  out- 
side the  state  in  order  to  solicit  their  comments  and  criticism.  This 
has  the  effect  of  calling  attention  to  our  program  and  gaining  the  bene- 
fits of  expert  criticism  which  will  improve  the  program.  A  distribution 
list  and  some  typical  comments  are  included  as  Exhibit  IV, 


6.0  EXTRA  REPRINTS  OF  MONTANA  SOLAR  DATA  MAOTJAL 

At  the  request  of  John  Omdorff ,  sixty  additional  copies  of  last 
year's  Montana  Solar  Data  Manual  were  reproduced.  The  initial  printing 
of  350  copies  had  been  exhausted.  This  activity  was  not  included  in  the 
original  contract  or  budget  but  sufficient  funds  were  gathered  from  var- 
ious budget  categories  to  pay  for  this  work.  These  copies  were  delivered 
to  the  Department  of  Natural  Resources  and  Conservation  in  June  198O. 


7.0  DATA  PROCESSING 

The  computer  hardware  and  software  were  altered  to  increase  the  ef- 
ficiency and  accuracy  of  the  solar  data  proceBsing.  We  now  have  a  hard 
copy  of  -the  data  for  each  solar  radiation  curve  as  it  is  read  into  the 
system.  This  change  of  format  gives  convenient  running  documentation  of 
the  data  base  which  is  kept  on  file.  A  sajnple  of  this  output  is  shown 
in  Exhibit  ¥. 


8.0  DATA  ACCURACY 

Our  Bolar  transducer  has  been  calitrated  on  two  separate  occasions 
by  Ed  Flowers  at  the  NCAA  laboratories  in  Boulder.  I  was  able  to  visit 
the  laboratory  and  discuss  the  results  during  October  1979«  The  NOAA 
calibration  sheet  is  shown  in  Exhibit  YI, 


9.0  CLEAR  DAY  MODEL 

An  improved  implementation  of  Douglas  Hoyt's  clear-day  radiation 
model  was  generated  during  the  past  year.  This  model  is  now  imple- 
mented on  the  same  computer  that  ws  use  to  compile  the  solar  radiation 
data*  In  previous  years  we  used  output  from  NOAA's  computer^  Now  we 
can  produce  original  comparisons  for  all  sites  right  in  the  office*  Oxnr 
version  has  been  compared  to  the  original  NOAA  version  for  Great  FallSj 
Montana  and  it  agrees  to  within  less  than  one  percent.  A  listing  of 
this  program  is  shewn  in  Exhibit  VII. 


10.0  FUIimS  FROJ-ECTS 

Ke  are  continually  looking  for  ways  to  improve  the  accuracy  and 
efficiency  of  our  task  of  measuring  and  documenting  the  solar  energy  re- 
source in  Montana*  We  are  currently  evaluating  the  feasibility  of  using 
a  solar  radiation  instrument  which  incorporates  digital  recording  in 
place  of  analog  chaJrt  recording*  We  hope  to  have  a  prototype  ready  for 
evaluation  this  yea^r. 

M®  have  initiated  a  comparison  of  our  solar  meastirements  with  two 
other  Montana  stations  that  are  part  of  the  National  Solar  Data  Program. 
Initial  comparisons  show  that  we  are  in  very  close  agreement.  This  com- 
psTison  will  be  extended  and  analyzed  during  the  coming  year. 


11»0  FINAL  BUDGET 

An  itemised  final  budget  report  is  shown  in  Exhibit  VIII.  The  ex- 
penditures at  the  end  of  the  grant  totaled  $34,794.00  leaving  a  net 
surplus  of  $971  that  was  not  used. 
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INTERPACINO  A  $10  CALCULATOR  FOH  DATA  ACQUISITIOS 


Charless  W.   Fowlkes,   Ph.D., 
POWLKIS   ENOINEERIHO 
31  Oardner  Park  Drive 
Bozeman,  MT     59715 


P.E. 


ABSTRACT 

Low  cost  instruments  that  can  be  used  to  make 
measurements  around  solar  projects  are  of 
interest  to  homeowners  and  solar  energy  ex- 
perimenters. For  solar  house  performance 
monitoring  it  is  important  to  know  the  total 
solar  radiation  intercepted  tiy  the  collector 
and  the  total  accumulated  degree  days.  A 
number  of  excellent  instruments  are  avail- 
able which  will  integrate  and  display  these 
quantities,  but  their  cost  may  exceed  the 
individual  experimenter's  budget. 

We  have  recently  built  a  circuit  which  re- 
quires less  than  $20  worth  of  parts  and  uses 
a  Texas  Instruments  hand  calculator,  costing 
about  $10,  for  a  read-out.  This  device  can 
be  connected  to  a  manufactured  or  home-made 
solar  radiometer  to  accumulate  and  display 
total  solar  radiation.  Another  application 
of  the  same  basic  circuit  uses  two  tempera- 
ture probes  to  effect  a  degree-day  integrator 
for  heat  loads  or  for  a  wood  otove  integrator. 
The  calculator  alone  can  be  used  as  an  event 
counter  to  display  pump  cycles,  door  openings, 
etc. 


1.   INTHODUCTIOH 

Some  hand  calculators  can  be  adapted  for  use 
as  a  low  frequency  counter.  Any  calculator 
which  will  increment  the  display  in  response 
to  a  single  key  stroke  can  be  used.  The 
Texas  Instruments  Tl-lOlO  used  in  our  experi- 
ments will  add  a  constant  to  the  display 
when  the  [=]  key  is  pressed.  For  example: 


Keystroke 

Display 

[+]  [1]  [=] 

I 

[-] 

■2 

[»] 

3 

or            "*''• 

C+]  [.25]  [-] 

.25 

[«] 

.50 

[-] 

.75 

etc. 

Pressing  a  key  on  a  calculator  closes  a  switch 
under  the  key.  The  procedure  used  to  inter- 
face the  calculator  as  a  counter  is  basically! 
(a)  Identify  the  key  which,  when  pressed,  adds 
a  constant  to  the  display  (the  [=]  key  in  the 
example  above),  (b)  locate  the  two  isides  of 
that  switch  in  the  calculator  keyboard  cir- 
cuitry, (c)  solder  a  wire  to  each  side  of  that 
switch  and  (d)  lead  the  wires  to  an  external 
switch.  This  external  switch,  which  is  in 
parallel  with  the  internal  keyboard  switch, 
will  now  increment  the  display,  effecting  a 
slow  speed  counter.  The  external  switch  can 
be  a  mechanical  switch  such  as  a  microswitch, 
a  relay  or  a  mercury  switch.  It  can  also  be 
an  electrical  switch  such  as  a  transistor, 
photo  transistor  or  Hall  effect  device. 

The  maximum  frequency  that  can  be  coimted  is 
determined  V  the  keyboard  scan  frequency  of 
the  particular  calculator,  which  is  usually 
between  5  and  100  Hz.  The  display  may  only 
respond  up  to  about  5  Hz  if  it  is  a  liquid 
crystal  type  (LCD) .  The  person  reading  the 
display  may  only  respond  up  to  one  to  two  Hz, 
which  sets  the  maximum  frequency  of  the  system. 

At  this  point  the  calculator  could  be  connect- 
ed to  a  switch  on  a  door,  a  vent,  a  movable 
insulation  system,  a  solar  circulating  pump 
or  a  ventilation  fan  and  it  would  display  or 
"count"  openings,  on-times,  etc.  In  this 
application  a  mechanical  counter  mi^t  be 
used  with  equal  cost  and  reliability. 

The  electrical  signals  from  temperature  trans- 
ducers, solar  cells,  pressure  transducers, 
heat  flow  sensors,  etc. ,  usually  consist  of  a 
voltage  or  current  which  varies  in  direct  pro- 
portion to  the  physical  quantity  being  sensed. 
There  exist  a  number  of  integrated  circuit 
chips  which  will  convert  a  do  input  voltage 
into  a  proportional  pulse  (or  switch  closure) 
frequency  on  the  output.  These  devices  are 
called  voltage-frequency  (V-P)  converters. 

By  adding  an  appropriate  V-F  converter,  we 
can  now  use  the  calculator  to  display  the 
integrated  total  of  a  variety  of  physical 
parameters.  The  steps  of  the  process  arei 
( a)  convert  the  physical  parameter  to  a  pro- 
portional voltage,  (b)  convert  this  voltage 
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into  a  proportional  frequency  of  switch  clo- 
sures and  (c)  count  the  switch  closurps  wit.^ 
the  modified  calculator  using  the  "constant 
add"  feature. 


2.      TRAJiSIXJCERS 

Solar  radiation  and  temperature  are  two   im- 
portant parametera  which  may  be  used  to   char- 
acterize the  performance  of  Bolar  heated 
tuildings.     A  low  cost  solar  radiation  trans- 
ducer can  be  constructed  with  a  silicon  cell 
and  a  wire-wound  shorting  resistor.     A  cell 
having  3.5  cm     area  and  a  1£L  shorting  resis- 
tor will  produce  about  80  mV  outpnit  in  bright 
Bun.     Properly  calibrated  and  tilted  toward 
the  sun  (at  typical  solar  collector  angles), 
this  device  will  have  an  average  accuracy 
approaching  1  5^  in  measuring  daily  total 
insolation.     Comnercial  solar  transducers  may 
also  be  used  with  this  calculator  circuit. 
Their  accuracy  will  be  somewhat  higher  but 
their  cost  will  be  much  higher. 

The  differential  temperature  circuit  included 
in  this  paper  uses  two  AI159O  transducers  manu- 
factured by  Analog  Devices.     These  transducers 
are  connected  to  produce  a  voltage  proportion- 
al to  the  difference  in  temperature.     Since 
the  calculator  will  onl;|r  count  or  integrate   in 
a  forward  direction,   one  probe  is  identified 
as  the  "hot"   probe  and  the  other  as  the   "cold" 
probe.     The  count  will  thus  be  proportional 
to  the  degree-hours  for  all  Thot>    Toold.      If 
the  hot  probe  is  placed  inside  the  house  and 
the  cold  probe  placed  outside,   the  calculator 
will   accumulate  actual  heating  degree-^ioura. 
If  the  probes  are  reversed,   the  calculator 
registers  cooling  degree-hours. 

3.     HARIWABE 

3.1     Interfacing  this  TI~1010 

Gently  piy  the  back  off  the  TI-1010.     Lay  the 
calculator  on   its   face  and  sli^tly  lift  the 


TABLE  1 
PARTS  LIST 

Hi,  R2,  ICK  1%  in  .22ufd  Mylar 

H4,  H9 

R3,  E7  lOOK  15^        C2  .0047ufd  Mylar 

R5  2.2K  1^        C3  .Olufd  Mylar 

E5,  Eg.  1.2UJ  15^       C4,  C5     lOufd 

Electrolytic 

RlOi  P-,1         28. 7K  1%       ICl  LM331 

RI2  2CSC  ten-       IC2  LMIO8A 

turn 

Dl  1M4148  IC3  ICL7C60 

VRl  m7805 


bottom  of  the  flexible  printed  circuit  to 
expose  the  terminals  on  the  keyboard  (see 
Pig.   1).     Solder  a  wire  to  the  farthest  left 
terminal.     This  wire  should  be  connected  to 
pin  3  of  ICl.     Solder  another  wire  to  the 
second  terminal   from  the  right.     This  wire 
should  be  connected  to  ground.     A   small  hole 
may  be  drilled  in  the  back  of  the  calculator 
case  to  allow  the  wires  to   exit  and  the  case 
to  close. 


Fig.  I.     Backside  of  TI-1010 


3.2    Circuit  Operation 

The  voltage  to  frequency  circuit  is  shown  In 
Fig.   2.     An  operational  amplifier  (IC2)   is 
configured  as  an   integrator  using  C2  as  its 
feedback  capacitor.     The  charge  rate  of  C2  is 
controlled  ly  -V^,  the  voltage  across  the 
load  resistor  (El)    and  R6.     R7  to  R12  provide 
an  offset  voltage  for  IC2.     When  C2  is  charged 
to  the  threshold  set  by  the  voltage  divider 
consisting  of  Rl  and  R2,  the  frequency  to 
voltage  converter  (ICI)   produces  a  switch 
closure  between  pin  3  and  ground.     A  TI-1010 
calculator  can  be  connected  to  count  and 
display  the  number  of  pulses.     R3  and  CI  are 
the  reset  timing  circuit. 

Pig.   2  also  shows  the  power  supply  for  this 
circuit.     A  power  supply  of  6f  do  or  greater 
is  required.     VRl  regulates  this  to  5?  do, 
and  IC3  converts  the  5?  dc  to  -JiT  dc.     Power 
consumption  of  this  circuit  is  only. 9. 2  ma 
(9-7  ma  with  two  temperature  transducers)   and 
can  be  powered  for  over  30  hours  of  inter- 
mittent use  with  a  SV  transistor  batteiy.     For 
long  tenn  use,   get  an  inexpensive  power  supply 
to  replace  the  battery.     A  6V  dc  supply  rated 
©  100  ma  (the  cube-plug  type  sold  for  transis- 
tor radios)   is  more  than  sufficient.     A  parts 
list  is  shown  in  Table  1. 

fig.  3  shows  a  completed  instrument  which  in- 
corporates both  the  silicon  cell  and  differ- 
ential temperature  inputs.     The  desired  input 
is  selected  with  a  switch.     (This  unit  was 
built  as  a  test  prototype  and  is  not  entirely 
practical  because  each  input  will  have  a 
diffarent   zero-offset.)     The  photograph  shows 
a  9V  battery  powering  the  circuit  for  porta- 
bility.    The  cube-plug  shown   in  the  photo- 
graph can  also  be  used  to  power  the  circuit. 


♦V 


I ^AAA— f — AAAA V^v\ — f-5\ , 


OUTPUT  TO 
CALCUUTOR 


TRANSCUCER 
IKPOT  '  Rl 


nri" 


Pig.  e.  Volta^-fre!]uan<7  o«nTart«r  and  powar  Buppljr  eiroults. 


Fig.    3.     Conpleted  prototype 

It  1b  a  nominal  6v  unit  which  costD  a  faw 
dollare  at  a  local  store. 

The  circuit  can  easily  be  wired  on  a  piece 
of  perf  board  about  the  size  of  the  calcula- 
tor.    Ne  put  the  prototype  circuit  in  a 
Btandai'd  audio  cassette  case  and  attached 
the  calculator  to  the  lid  of  the  case.     This 
was  convenient  as  the  lid  of  the  case  could 
ba  opened  to  aooess  the  oiixiuit  for  oheokout 


Fig.   4.     Cane  opened  to  show  prototyp*  olroult. 

and  zero  adjuetment.     A  photograph  with  the 
cover  opened  la  shown  in  Pig.   4.     Whan  using 
a  solar  cell  as  an  input  device,  Rt   must  be 
approximately  . 3Jl  .     The  solar  cell  muat  be 
connected  with  the  positive  terminal  to 
ground.  Fig.   2. 

When  using  temperature  transducers  as  input 
devices,  Rj^  should  be  selected  as  followsi 
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For  temperature  differentiale  of  O-50  c, 

Rl~  %Jt  . 

For  tf  ■■^erature  differentials  of  0-25°  C, 
Rl-^  '-■  .Ti.  ■ 

The  "cold"  temperature  probe  must  have  +3?   dc 
on  the  positive  lead  and  the  "hot"  probe  must 
have  -5V   dc  on  the  negative   lead.     Fig.    5 
BhowB  thie  input  circuit. 


"cold" 


"hot" 


Pig.   5'     Differential  temperature  input 
circuit. 


4.     USIKO  THE  CIRCUIT 

normal  precautions  should  be  taken  before 
connecting  power  to  the  circuit.     After  the 
circuit  has  been  built  ajid  de-bugged,   the 
following  procedures  will   prepare  the  instru- 
ment for  use. 

4.1     Adjusting  the  Offset  Voltage 

With  silicon  solar  cell  as  input  devioei 


i] 


4) 


Short  the  input  load  resistor  R, , 
Connect  a  voltmeter  to  pin  6  of  IC2. 
The  charging  ramp  of  C2  should  be  ob- 
eerved.     The  voltage  should  be  between 
approximately  -.6  and  S.JIf  do. 
3)  If  the  voltage  is  constant,   adjust  H12 
until  the  ramp  begins  to  cycle.     Hots 
the  value  of  the  peak  voltage. 
Slowly  adjust  H12  imtil  the  ramp  volt- 
age holds  steady  at  about  O.IV  below 
the  peak     (i.e.   if  the  peak  is  2,5»f, 
adjust  R12  irntil  the  voltage  is  steady 
at  2.4V. 

With  AD59O  temperature  probes  as  input  device: 

1)  Ismierse  both  of  the  AD590  temperature 
probes  in  an  ice  bath  e^d  allow  their 
temperature  to  stabiliKS. 

2)  to  4)  This  procedure  the  same  as  above. 

4.2    Calibration  with  Silicon  Solar  Cell 

If  you  do  not  have  a  refiBremce  lyranooetar, 
the  following  clear  day  procedure  should 
allow  calibration  to  within  1  I05J  or  better. 


Calibration  should  be  done  on  a  clear  day 
within  one  hour  of  solar  noon. 

1)  Set  the  offset  voltage  as  described  in 
Section  4.1. 

2)  Turn  on  calculator  and  key  in  [+]  [1]. 

3)  Hold  the  Bolar  cell  normal  to  the  sun 
(directly  facing  the  mm). 

4)  Record  the  calculator  count  and  begin  a 
one  minute  test.     After  one  minute 
record  the  new  count.     IJOTEi     If  the 
calculator  count  is  too  rapid  to  read, 
the  size  of  the  load  resistor  (R,) 
should  be  decreased, 

5)  Look  up  the  Beam  Radiation  at  solar 
noon,  I,   for  your  latitude  from  Table  2. 

6)  Add  1%  to  I       for  each  1000  feet  of 
elevation  above  saa  level.     Add  3J5  for 
a  very  clear  (arid  or  mountain)   atmos- 
phere.    Subtract  %  for  coastal  or 
polluted  atmospheres. 

7)  The  solar  calibration  constant  is 


(l))iH/ni^ 
count/rain 


hr 

60  min 


lOOOH 


"      1^.7  3c  I  Wh-m-E 
£i  count 

Ezamplei     In  one  minute  50  counts  were 
recorded.     From  Table  2  for  the  month 
of  June  and  40°  North  Latitude,   1-0.88. 
Elevation  is  3000  feet,  therefore: 

I  -  0.88  +  (0.03  I  0.88)   .  0.906 

An  additional  ^  is  added  for  a  very 
clear  atmosphere! 

I  -  0.906  +  (0.03  2  0.S8)  . 

0.931  kH-oi-2 
Constant  >  16. 7  x  O.93I  Wh-m"^ 
60  counts 
-  0.26  Vh-m-^-ct-l 
8)  To  use  the  device,  ke(y  in  [+]  [.26]. 
The  display  will  now  accuimilate  Wh-o-8 
of  reudiation  falling  oa  the  silicon 
cell.     Any  oonvsnient  units  can  be  used 


TABLE  2 
HORKAL  BEAM  RADIATICW  AT  SOLAR  BOOM  (BI-o-2) 

24**  UT      32°LAT-     AO°ht^      48°  LAT 


JAN 

1.01 

0.98 

0.93 

0.64 

FKB 

1.02 

1.00 

0.97 

0.93 

KEftR 

1.00 

0.99 

0.97 

0.94 

i?R 

0.94 

0.94 

0.92 

0.91 

KAY 

0.91 

0,90 

0.90 

0.88 

JUH 

0.89 

0.88 

0.88 

0.87 

JUL 

0.88 

0.88 

0.87 

0.86 

ADO 

O-90 

0.90 

0.88 

0.86 

SEP 

0.95 

0.93 

0.91 

0.86 

OCT 

0.98 

0.96 

0.93 

0.88 

HOV 

0.99 

0.96 

0.91 

0.82 

lEiC 

1.00 

0.96 

0.96 

0.79 
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V  simply  keying  in  the  appropriate 
conBtant. 

4«3  Calibration  with  Temperature  Probee 

1)  Set  the  offset  voltage  ae  described  in 
Section  4.I. 

2)  Turn  the  calculator  on  and  key  in 
[+],  [1]. 

3)  Select  two  water  bathe  with  tempera- 
tures differing  approximateljr  'i^  the 
amount  corresponding  to  the  load 
resistor  (Ri)  being  used  (i.e.  25°  C 


4) 


for  R, 


IK 


Immerse  the  "cold"  probe  in  the  cooler 
bath  and  the  "hot"  probe  in  the  warmer 
bath  and  allow  to  stabilize. 

5)  Eecord  the  calculator  count  and  begin  a 
one  minute  test.  After  one  minute 
record  the  new  count.  NOTKj  If  the 
calculator  count  is  too  rapid  to  read, 
the  temperature  differential  or  the 
size  of  the  load  resistor  (R^)  should 
be  decreased.  If  the  count  appears 
very  sluggish  (i.e.  only  10  counts  per 
minute) ,  the  temperature  baths  may  have 
too  small  a  temperature  differential  or 
Rj^  should  be  increased. 

6)  Record  the  temperature  of  the  "cold" 
bath  (T(.)  and  the  temperature  of  the 
"hot"  bath  (Th). 

7)  The  calibration  constant  =   Tjj  -  Tq 

oount/min 

Example;     In  one  minute  60  counts 
were  recorded.     The  tranperature  of  the 
cold  bath  (Tq)  was  0°  C   and  the  temp- 
erature of  the  hot  bath  (Tji)  was  25°  C. 

Constant  ^     (25-0)°C         =0.4l67°C-min/ot 
60  c6unts/min 

8)  To  use  the  device,   clear  the  display 
and  key  in  [+]  [.4I67].     The  display 
will  now  accumulate  °C-min  when  Thi^  Tg. 
The  constant  can  be  converted  to  display 
other  convenient  units.     For  example, 
0.4167  X  1.8/60  =.  0.0125  °P-fer/oount. 


5.     RECOmmPATIOHS 

He  have  built  two  of  these  circuits  and  they 
have  worked.     The  design  should  be  considered 
as  a  prototype  and  the  user  should  be  pre- 
pared for  possible  bugs  that  we  may  have 
overlooked.     The  Tl-lOlO  calculator  that  we 
used  for  the  prototypes  will  autoraatioally 
turn  itself  off  after  several  minutes  of  no 
keyboard  entiy.     Another  calculator  without 
this  feature  would  have  to  be  used  for  long 
time  measurements. 


anyone  wanting  to  explore  other  applloatloDB 
which  use  calculators  to  display  physical 
measurements.     Reference  (3)   disousses 
details  of  using  low-cost  silicon  cells  to 
measure  insolation.     References  (4)   and  (5) 
discuss  performance  monitoring  of  solar 
houses  using  inatrumants  of  moderate  cost  and 
may  be  of  interest. 

The  parts  for  this  circuit  can  be  ordered 
from  Jameco,  TriTek,  Electrolabs,   or  other 
selected  mall  order  electronics  parts  houses. 
The  addresBes  of  electronic  parte  suppliers 
can  be  found  in  advertisements  in  most  hobby 
electronics  or  computer  magazines. 
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ABSTRACT 

Per  tho  past  four  yearB  in  Montana  we  hawa 
operatad  a  network  of  30  low-cost  solar 
radiation  measui-ing  laiitruments   (1) .    At 
this  point  the  data  liosie  consists  of  over 
40,000  data  days.     Tho  hasio  elemsnt  of  data 
Btorad  on  the  computer  is  the  half-hour 
average,  radiant  intensity  on  a  60°  tiltad 
surface.     Tba  sensors  are  silicon  calls 
mounted  on  a  flat  plate  and  connected  to  a 
strip  chart  recorder.     The  cost  of  aach 
instrument  at  the  'begiiiming:  of  the  SUSK 
prograan  was  ahout  $200. 


1.      IWRoajCTIOM 

The  instruments  are  mounted  on  public  high 
schools  throu^iout  Montana.     This  program 
uses  volunteers,  typically  high  school 
teachers,  who  help  maintain  the  instruments 
in  the  field.     The  strip  chart  records  are 
mailed  hy  'the  volunteers  to  our  office 
where  they  are  traced  on  a  digitiaer  for 
input  into  a  microcomputer. 

The  output  of  our  solar  meamipemant  program 
is  publlshad  annually  in  tha  Montana  Solar 


9    Continuous  Data  Recording  Station 
O  Haziual  Data  Station 

Pig.   1.     Tha  SIMM  Network  (Solar  Insolation 

MeaBurement  IJontana) 

Maintaining  the  aoouraqy  and  continuity  of 
93tperimental  data  is  a  difficult  task;  due 
to  tha  large  numher  and  state-wide  distribu- 
tion of  tha  sites    (Fig.    1).   To  visit  all   30 
instrumented  sites  requires  travelling  over 
3,000  miles.     Since  the  instruments  are 
home-made,  their  performance  ntust  be  care- 
fully and  scientifically  documented  to  re- 
late the  radiation  data  to  accepted  calibra- 
tion standards.     Tho  content  of  this  p^er 
relates  our  approach  to  these  tasks.     Spaoi- 
fioally,   tha  results  of  (a)  NOAA  calibrations, 
(b^   aida-ty-aide  pyranometer  comparisons, 
(o)  clear  day  theoretical  comparisons  and  (d) 
oomparison  with  independent  measurements  from 
the  National  Solar  Data  Hotwork  are  discussed. 
The  results  of  each  appiroach  in  asaeoaing 
accuracy  and  maintaining  data  integrity  of 
the  3IHU  instruraants  is  discussed. 


montaoa 


solar    data 


manual 


Fig.  2.  Cover  of  Montana  Solar  Data  Manual, 
illustrating  insolation  records  ( including 
eclipse  on  2/26/79) 
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Data  Manual    (2) .     TliiB  manual  addressas  an 
av3'-:ice  that  ranges  from  technical  and 
J"       •  sional  people  to  individual  hoaie- 
c./i^xs,  backyard  inventors  and  solar  anthu- 
siaata.     The  manual  reporta  insolation  data 
and  cloudy  day  statistics,     Climatological 
data,   solar  design  methods  and  a  selected 
bibliography  are  also  included  in  tlie  majnual. 
TSiie  manual  is  distributed  f^e  by  the  spon- 
soring agancy.  Alternative  Renewable  Energy 
Sources  Prograa,   Montima  Department  of 
Natural  Resources  and  Conservation.     About 
500  copies  of  the  manual  have  been  distribu- 
ted to  libraries,   universities  and  designers. 
Another  500  copies  of  a  revised  and  updated 
manual  are  being  prepared  for  distribution 
during  the  summer  of  1981.     The  SIMM  network 
is  in  its   fifth  year. 

2.     BATA  REPOCTICM  KRSOR 

S'trip  chart  d&ta  is  processed  117  tracing  the 
curves  on  a  digitizer  platan  connected  to  a 
microcomputer.     The  resolution  of  the  platen 
is  0.1  nm  and  each  curve  is  represented  Ijy 
200  to  400  coordinates.     The  reading  error  due 
to  resolution  is  lass  than  1%,     The  primjuy 
source  of  data  reduction  error  is  the  repeat- 
ability of  the  ourve  tracing  itself.     Repro- 
ducibility teste  have  been  run  on  over  1,000 
data  days.     The  reproducibility  of  the  month- 
ly average  insolation  is  always  ±  25S  imd 
often  within  1^. 


3.     TRASSPUCER  ERRORS 

The  SIMH  silicon  cell  transducer  is  covered 
by  a  glass  plats  potted  in  silicone  sealant. 
Tilt  calibration  tests  have  shown  that  the 
output  of  the  sum  transducer  falls  off  at 
large  angles  of  incidence.     This  is  due  pri- 
marily to  the  glass  cover  plate  and  is 
normally  tanned  "cosine  error". 

The  measured  cosine  response  of  the  SIHM 
transducer  was  input  into  a  computer  simula- 
tion modal  to  predict  the  effect  of  this 
error  on  the  daily  total  insolation  thrtiu^- 
out  the  year    (Table  1) .  Ezperimental  data 
presented  later  in  this  report  indicates 
that  the  suimer  errors  are  less  than  those 
predicted  ly  the  modal. 


4.     CIZAR  BAT  METHOD 

All  SIKM  tranitducers  are  calibrated  before 
installation.     The  reference  instrunects  have 
included  Kipp  pyranoaeters,  Li-Cor  radio- 
meters and  a  NCAA  calibrated  SIKH  transducer. 
To  allow  for  drift  aai  aging,  which  may  occur 
while  the  instrument  is  on  site,  we  malce  use 
of  a  clear  day  correction  technique  similar 
to  the  methods  used  to  correct  the  SOUJET 
data  bass. 

The  essence  of  our  clear  day  prooedore  is  to 
compare  the  measured  noon  fltix  to  a  value 
predicted  t|y  a  theoretical  clear  day  model. 
We  have  u£>8d  the  ASHRAS  model  and  a  nodsl 
pwre-ided  Tecf  Douglas  Hoyt  of  WOAA     {  3) .  As 
part  of  our  data  reduction  procedure  we  pick 
off  the  noon  maximum  on  all  clear  days  and 
compute  a  theoretical  value  for  that  eama 
day.     If  thare  is  a  consistent  difference  of 
more  than  n  few  percent,  the  ezperjjneotal 
data  for  that  period  is  all  scaled  to  match 
the  theoretical  value.     Figure  3.Bhows  a 
graph  of  experimental  and  theoretical  clear 
day  values  for  two  stations  for  1978. 


120  160     200  240     280  320   360 


Pig. 3a.      Great  Palls  data  for  I978 
showing  two  clear-day  models  with  data 


1.1  I- 


.  scale  factor  '*--a^<Br''^ 


0     40     80       120  160     200  240     280  320  360 

Fig.3b.GlondivB  data  for  1978 
'illustrating  use  of  soale  factor 


TABLE  1. 
CALCULATED  EFTECT  OP  COSINE  ERROR    (») 


MOUTH 

CLEAR 

iwr 

AVERAGE  BAT 

use 

-2 

-2 

JAK       NOV 

-3 

-£ 

FEB       OCT 

-4 

-3 

KSB       SEP 

-4 

-3 

APR       AUG 

-6 

-4 

HAT       JDL 

-8 

-5 

jhh 

-8 

-5 

Experimental  clear  day  maximums  from  an 
Eppl^  PSP   (used  in  the  National  Solar  Data 
ITotworfc)   are  compared  to  the  theoretical  Hoyt 
values  in  Table  2.     The  measured  maximums 
fluctuate   sbout  the  theoretical  maximum. 
This  coul!'  be  due  to  fluctuations  in  the 
atmosphere  and  errors  due  to  temperature 
effects,   calibration  and  dirt  on  the  instru- 
ment.    The  absolute  value  of  this  fluctua- 
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■tton  for  i960  in  Billines ,  MT  ayerage  2.3%. 
SIKH  data  are  -typioall^  matched  to  the  Hoyt 
model  to  within  3^  anaual  average  deviation 
(1,  2).  This  data  thiui  mipporta  the  clear 
day  method  used  in  the  SIMS  network, 

TABLE  2., 
CLEAR  MI,   NOON  WAX,   (kW-m  ),  BILLIHCS  I98O 


a  dozen  points  that  are  scattered.  These 
disorepanaies  probably  occurred  on  days  when 


DAT 


HOTT 
SOAA 


NSM 


POWLKES 

(sum) 


6 

19 

20 
28 
29 
54 
55 
36 

95 
97 
104 
105 
107 
108 

iir 
117 
118 
127 

128 
142 
143 
145 
161 
165 
174 
177 
178 
179 
180 


0.92 
0.95 
0.95 
0.98 
0.98 
1.06 
1.06 
1.08 
1.05 
1.05 
1.02 
1.02 
1.01 


1 

1 
0, 
0, 
0, 
0, 
0, 
0, 
0, 
0.86 

0,85 
0.85 
0,85 
0.85 
0.85 
0,85 


,01 
.00 

,98 
97 
94 
94 
89 
89 
88 


0.98 
0.96 
1,01 

1.07 

1,17 
1.04 
1.07 
1.07 
1.03 
1.05 
1.02 
0.97 
1.00 
1.00 

0.96 
0.95 
0.93 
0.90 
0.92 
0.84 

0.86 
0.86 
0.84 
0.83 
0.82 
0.84 
0.84 
0,86 
0.86 


0.98 
1.00 
1.01 

1.01 
1.00 
1,01 

1,07 
1.05 
1.03 

1.00 
1.00 

0,98 


01 

99 
99 
98 
95 
92 
94 


0,90 
0,90 

0.92 
0.83 
0.85 
0.83 
0.83 
0.84 
0.85 
0.87 


5-     NATICMAL  SOLAR  DATA  HBWORK  COMPARISON 

The  Billings  Shipping  solar  system  in  Billings, 
Montana  was  instrumented  three  years  ago  for 
performaaoe  monitoring  ao  part  of  The  National 
Solar  Data  Network  (NSns) .     alobal  insolation 
is  measured  hare  by  an  Eppiey  PSP  pyranometer 
tilted  at  50°.     Daily  totals  and  nmximma  are 
reported  each  month  in  Environmental  Data  for 
Sites  in  the  National  SoTar  Data  Network     (A). 

During  I98O  there  were  337  days  of  siJiiultaa- 
eous  data  from  both  the  NSDI   (Eppiey  PSP)  and 
the  SIMH  transducer  at  Billings  Hi^  School. 
The  NSISf  50°  tilt  data  was  converted  to  a  60° 
tilt  using  the  ratio  of  ASHRAE  clear  day 
totals.     This  simplified  method  would  tend  to 
over-correct  the  KSW  data  by  a  few  percent 
in  the  sunmer.     Since  the  correction  itself 
averages  only  about  t  3?Sf  the  error  intro- 
duced is  not  large     (5)  . 

A  scatter  plot  of  daily  total  insolation  for 
SIMM  vs.  NSai  for  I98O  is  shown  in  Figure  4. 
For  337  common  data  days,  the  linear  oorrala- 
tion  coefficient  was  0,990.     There  are  about 


337  Days 
1980 
Billings,  HT 


Linear 
correlation 
coefficient 
0.990 


2  4  6  8 

Eppiey  PSP 

Pig.   4.   Scatter  plot,   daily  totals,   Whra^ 


TABLE   3. 
DAILY  TOTALS,    (M-J-ra^  )  ^  BILLINCS,  NOV.    I98O 


NSIN 

NSKf 

POVJLKES 

LESS 

DAT 

(PSP) 

SIMM 

SIMM 

1 

2.59 

2.66 

-.07 

i 

2.18 

2.13 

.05 

3 

5,70 

5.72 

-.02 

4 

1.92 

1.59 

.33 

5 

6.07 

6.16 

-.09 

6 

1.51 

1.25 

.26 

7 

1,06 

1.00 

.06 

a 

4.41 

4.78 

-.37 

9 

1.06 

1.06 

.00 

10 

2.42 

2.47 

-.05 

11 

1.23 

1,19 

.04 

12 

0.65 

0.63 

.02 

13 

0.87 

0.47 

.40 

14 

5.27 

5.06 

.21 

15 

5.80 

5.91 

-.11 

16 

5.74 

5.84 

-.11 

17 

5,87 

5.81 

.06 

18 

2.92 

2.84 

.08 

19 

1.65 

1.41 

.25 

20 

3.76 

3.84 

-.08 

21 

1.82 

1.84 

-.03 

22 

0.58 

0.50 

.08 

23 

5.75 

5.63 

.13 

24 

2.28 

2.22 

.06 

25 

0.85 

0.53 

.31 

26 

5.17 

4.63 

•  55 

27 

0.94 

0.81 

.13 

28 

5.80 

5.31 

.48 

29 

1.82 

1.72 

.11 

30 

0.39 

0.38 

.02 

AVG 


2.94 


2.85 


3.07fo 
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enow  or  trou't  covered  one  •transducer  part  of 
that  day.     Nslther  Instrument  ie  regalarljr 
cleaned  or  cleared  of  snow.     The  remaining 
scatter  is  primarily  (^  n  to  differences  in 
apectr&l  response,   ca     sration  error,   oosine 
response  and  teoipersture  coefficients  (annual 
temperature  extraoes  span  over  70'  C). 

The  daily  totals  for  the  month  of  Hovenber, 
i960  are  shown  in  Tahle  3«     The  daily  average 
insolation  for  Novemher  for  the  SIMM  instru- 
Bont  is  3.07^  less  than  the  NSXR  instroment. 
This  difference  closely  oorrcapondjo  to  the 
-3j5  error  predicted  for  this  wnsth  due  to 
the  low  oosine  response  of  the  SDH  trans- 
ducer. 

The  avara^  dally  totals  for  aach  inatruosnt 
are  shotm  for  each  eonth  in  1^60  in  Tahlc  4< 
The  SmK  Instrument  is  normally  louar  ^  a 
few  percent,  which  is  prohahly  due  to  cosine 
error.     (The  Decsmber  a^era^  is  faaltad 
hecause  the  SBSH  instrument  was  not  rsoordisg 
data  during  a  week  with  several  clear  days.) 


POtfliKES 

SBW 

3.29 
3.29 
4.01 
5. 18 
4.74 
4.82 
5.15 
4.73 
4.72 
4.33 
2.85 
1-88 


TABIi  4. 

MOTTHLI 

AVERAOE  IDfTT.T  TOTAL  (W 

HSOf                      I 

MOTTS 

(PS?) 

JUS 

3.13 

FEB 

3.29 

BUR 

4.14 

APR 

5.23 

MAT 

4.77 

JU» 

4.84 

JUL 

5.25 

ADO 

4.75 

SEP 

4.80 

OCT 

•4.38 

NOV 

2.94 

TEC 

2.05                  « 

'^missing  8  to  10  days 


6.     SHE  BT  SIDE  CCMPARI30N3 

Side  tiy  side  compariaons  were  run  for  ei^t 
months  during  the  winter  of  I98O  and  sprijig 
of  1981  to  compare  (a)   a  Kipp  Hyranometer, 
(h)   a  Li-4or  200S   silicon  cell  radiometer, 
and  (c)  the  SIKH  silicon  cell  radiometer. 
All  transducers  were  mounted  at  0°  azimuth 
and  60°  tilt  and  connected  to  a  microcorapater 
'b'uied  data  acquisition  systea  (SAS).     The  BAS 
:.■':.  \  resolution  of  12  hits,   on  accuracy  of 
L^Uur  than  t.  0,%;   each  transducer  and  -the 
ambient  air  temperature  were  sampled  once 
each  second.     These  readings  were  averaged 
over  10  minute  intervals  and  the  average  data 
stored  on  magnetic  tape.     This  data  base 
consists  of  13  X  10°  side  by  aide  readings  of 
each  transducer. 


Table  5  shows  a  typical  portion  of  this  data 
base  condeneed  into  daily  total  insolation. 
The  Kipp  pyranomstar  was  recalibrated  at  the 
KOAA   facility  in  1979  and  this  calibration 
factor  was  uBud.    The  faotoiy  calibration 
supplied  for  the  Li.<;or  was  used  and  the  SIXH 
transducer  was  calibrated  to  the  Kipp  on  a 
clear  ds^  near  solar  soon  acoording  to  onr 
standard  procedure.    Ho  correction  was  made 
for  ambient  temperature,  which  varied  betweeo 
-18°  C   and  ♦IT"  C  during  this  period,     Th® 
monthly  totals  of  the  Kipp  and  the  Li-^or 
agree  to  within  t  1%  while  the  SIHK  tnms- 
duocr  is  39(  below  the  Kipp.    This  lower  read- 
ing agrees  with  the  predictions  of  the  oosine 
error  model  discussed  previously  (Table  I), 

TABLE  5. 
BOZBUH  BAILY  TOTALS*  (iWh-o^) 

lorajcBiE,  1980 


DAI 


KIPP 


LI -COR 


aim 


I 

2.80 

2 

1.95 

3 

1.01 

4  ' 

2.04 

7 

0.42 

a 

4.21 

9 

1.39 

10 

1.40 

li 

0.54 

12 

1.14 

13 

5.67 

14 

*5.33 

15 

6.66 

16 

6.62 

17 

4.13 

18 

1.10 

19 

0.93 

20 

3.83 

21 

3.04 

22 

0.73 

23 

5.43 

24 

3.71 

25 

0.71 

26 

4.44 

27 

1.74 

29 

2.01 

30 

1.54 

A.TO 

2.66 

2.64 

2.72 

1.94 

1.88 

0.99 

0.93 

2.07 

1.96 

0.44 

0.39 

4.22 

4-13 

1.39 

1.31 

1.32 

1.30 

0.55 

0.53 

1.15 

1.12 

5.62 

5.44 

6.60 

6.36 

6.55 

6.38 

6.61 

6.43 

4-40 

4.24 

1.03 

1.05 

0.89 

0.86 

3.77 

3.68 

3.02 

3.00 

0.68 

0.66 

5-40 

5.28 

3.61 

3.54 

0.64 

0.63 

4.43 

4.28 

1.66 

1.66 

1.94 

1.93 

1.45 

1.44 

2.64 


2.57 


The  instantaneouis  readings  of  the  three 
transducers  during  mid-day  in  clear  weather 
were  always  within  a  ^  ^  band  and  often 
were  within  ±  IfL     The  Kipp  and  the  Li-Cor 
were  typically  within  i  0.01  l*-«j^     for  the 
entire  day,  while  the  SIMM  transducer's  out- 
put would  typically  be  0.05  W-ai^     low  in 
early  morning  and  late  afterc  r>n  (the  oosine 
error). 
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7.     Snai  HOT/,'ORK  RSLIABILm 

Each  summer  a  teolmician  visits  all   30  sta- 
tions to  check  the  transducer  mounting  and 
the  condition  of  the  wiring  and  the  strip 
chart  recorder.     If  the  system  breaks  down 
during  the  yeaj  the  recorder  and/or  trans- 
ducer sire  sent  to  us  hjr  the  volunteer.     V'e 
repair  or  replace  the  unit,  test  it,   and 
send  it  back. 

During  the  period  hetviresn  January,   1977  and 
December,   I980,  there  were  38,956  possible 
station  data-days.     Vn  recorded  32,425  data- 
days,   or  83^  of  the  total.     The  IT^S  lost  data 
was  primarily  due  to  recorders  running  out  of 
paper  (mostly  during  the  sunmer  months). 
During  four  years  of  operation  two  recorders 
were  stolen,  two  transducers  were  vandalized, 
six  motors  burned  out  and  there  were  about 
two  dozen  instances  of  miscellaneous  mechani- 
oal  or  electrical  failure.     Data  was  lost 
several  times  when  wires  were  cut  or  power 
was  turned  off  during  remodeling  worfc  at  the 
school.     Some  data  has  been  lost  in  the  mail 
and  in  two  instances  strip  charts  were 
destroyed  by  a  janitor's  floor  polisher.     Some 
suspicious  data  has  been  rejected  after 
arriving  at  the  office. 


8.     COHCLUSIOKS 

•  The  clear  day  correction  method  provides 
an  economical  way  to  maintain  the  calibration 
of  the  network.     Currant  comparison  data 
indicated  that  this  correction  puts  the  data 
into  a  i  5^  accuracy  band. 

«       The  data  reooveiy  of  83^  would  be  diffi- 
cult to  improve  in  this  volunteer  opsrated 
netvioric.     Missing  data  is  not  too  important 
due  to  the  high  density  of  the  30  meaouring 
stations. 

•  Extensive  compariaon  data  shows  that 
errors  due  to  spectral  and  temperature  effects 
are  well  within  the  ±  %  error  band.     The  most 
serious  fault  with  the  SUM  trsasduoer  is  the 
poor  cosine  response  which  produoes  a  aystam- 
atio  error  of  about  -3=^  during  the  heatiag 
season. 
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«       This  data  base  has  been  made  available 
to  the  public  along  with  supplementary 
material  in  an  effort  to  improve  the  quality 
of  solar  design  in  Montana- 

e  Measuring  insolation  on-  a  6O0  tilted 
surface  reduces  the  net  error  in  the  data 
when  used  for  solar  desigp. 
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MEASUREMENTS  OF  SOUR  RADIATION   IN  THE  UPPER  HEMISPHERE 
USING  AN  ARTICULATED  RADIOMETER 


CharlessW.  Powlkes,  Ph.D.,  P.E. 
Fowlkes  Engineering 
31  Gaxdner  Park  Drive 
Bozeman,   Montana     59715 


ABSTRACT 

The  goal  of  this  project  is  to  measure  solar 
radiation  (quasi)  simultaneoUBly  at  144 
positions  of  tilt  and  azimuth  in  the  upper 
hemisphere.  The  apparatus  consiBts  of  a 
solar  radiometer  fixed  to  an  articulating 
mechanism  driven  by  stepping  motors  which 
are  controlled  hy  a  low-cost  microcomputer. 
At  each  15°  increment  of  altitude  and  azimuth, 
the  radiant  flux  is  sensed  by  the  radiometer 
and  the  resulting  electrical  output  digitized 
and  stored  in  the  computer.  The  upper  hemi- 
sphere is  scanned  once  each  30  seconds.  At 
the  and  of  each  hour,  the  solar  radiation 
measured  in  each  direction  is  averaged  and 
the  result  stored  on  magnetic  tape.  The  data 
base  consists  of  hourly  averages  of  the  solai" 
flux  on  each  of  144  planes.  This  paper 
describes  the  experimental  apparatus  and  pre- 
sents preliminary  data. 


1.0  INTRODUCTION 

The  majority  of  available  solar  radiation 
data  consists  of  global  radiation  measured  on 
a  horizontal  plane.  A  horizontal  orientation 
is  normally  selected  by  meteorologists  because 
they  are  primarily  interested  in  a  "horizon- 
tal" solar  collector,  the  earth.  A  solar  de- 
signer, however,  will  normally  tilt  the 
collector  plane  between  40  and  90°  up  from 
the  horizontal  to  increaise  the  average  radi- 
ant flux  density  intercepted  by  the  collec- 
tor. To  meet  the  needs  of  solar  designers, 
a  number  of  empirical  rules  have  devised 
relating  solar  flux  on  a  horizontal  surface 
to  the  solar  flux  on  a  tilted  surface.  A 
complete  review  of  the  literature  of  "tilt 
correlation"  studies  is  beyond  the  scope  of 
this  paper,  but  samples  of  important  work  in 
this  area  are  included  in  References  (l),  (2), 
(3),  and  (4). 

A  comparison  of  several  popular  tilt  correla- 
tion algorithms  was  initiated  at  the  Solar 
Energy  Research  Institute  (SERI)  in  1978. 
This  comparison  revealed  discrepancies  among 
the  various  algorithms  which  stimulated  the 
initiation  of  experimental  measurements  of 


insolation  on  tilted  surfaces  (5).  Insola- 
tion data  has  subsequently  been  collected 
at  the  SERI  field  site  from  an  apparatus 
which  consists  of  a  series  of  pyranometers 
facing  south  and  tilted  at  angles  between 
90  and  0"  at  10°  increments.  This  array  of 
measuring  instruments  included  Eppley  PSP 
pyranometers  and  Li-Cor  IJ:-200S  silicon  cell 
radiometers  mounted  side  by  side.  Data  from 
these  experiments  will  probably  be  published 
during  the  coming  year.  Data  from  siraultane- 
OUB  meaaurements  on  six  planes  is  given  in 
Reference  (7). 

Silicon  cells  have  a  non-uniform  spectral 
response  which  will  produce  errors  in  the 
indicated  radiant  flux  in  some  situations  (6). 
Preliminary  data  from  the  SERI  experiment 
showed  close  agreement  between  the  global 
solar  radiation  indicated  by  the  pyranometers 
and  the  solar  radiation  indicated  by  the 
silicon  cell  based  radiometers  (5),  (3).  Our 
data  from  extensive  side— by-side  pyranometer 
and  silicon  cell  radiometer  tests  have  also 
shown  good  overall  agreement  (9).  Since  the 
silicon  cell  radiometer  exhibits  both  ade- 
quate accuracy  and  time  response,  it  appeared 
feasible  to  use  this  transducer  in  the  scan- 
ning mode  to  gather  insolation  data  (quasi) 
simultaneously  at  a  large  number  of  angles  of 
incidence. 


'2.0  SCANtJER  APPARATUS 

Figure  1  shows  a  schematic  of  the  overall 
apparatus.  A  silicon  cell  radiometer  is 
attached  to  a  stepping  motor  which  will 
adjust  its  altitude.  This  assembly,  in  turn, 
is  mounted  on  a  second  stepping  motor  with  a 
vertical  axis  which  can  adjust  the  azirmith  of 
the  assembly.  Both  stepping  motors  are 
connected  to  appropriate  drive  electronics 
which  are,  in  turn,  interfaced  to  a  micro- 
computer. The  altitude  and  azimuth  of  the 
solar  radiometer  can  now  be  adjusted  in 
response  to  commands  issued  by  the  micro- 
computer. The  output  of  the  solar  radiometer 
is  connected  to  an  a/D  converter  (Figure  1) 
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Fig.  1.  Schematic  diagram  of  computer  control  and  interface. 


which  is  interfaced  to  the  microcomputer. 
This  experimental  arrangement  allows  the 
computer  to  move  the  sensor  to  a  specified 
azimuth  and  altitude  and  to  read  the  corre- 
sponding solar  radiation. 

The  solar  radiometer  chosen  for  this  experi- 
ment was  a  Li-Cor  LI— 200S.  The  mechanical 
arrangement  is  shown  in  greater  detail  in 
the  diagram  in  Figure  2,   In  this  diagram, 


Altitude 
motor 


Fig.    2.     Details  of  actuator. 

the  radiometer  is  mounted  to  the  shaft   of  the 
altitude   stepping  motor.     This  assemlily,    in 
turn,    is  mounted  to  the  shaft  of  the   azimuth 
stepping  motor.     The   altitude  stepping  motor 
can  rotate  through   180°   in  7.5°  increments 
(from  a  vertical  tilt   through  horizontal  to   a 
vertical   tilt   in  the  opposite  direction) , 
Figure  2.     The  azimuth  motor  can  rotate 
through    >180°   in   15°   increments. 


3.0     COMPUI'ER   INTERFACE 

The  microcomputer  used  with  this   apparatus   is 
a  Radio  Shack  TRS-80  with   16K  memory   and 
Level   II  BASIC.      This  machine  was   chosen  be- 
cause we  have  used  several  TRS-80 's   for  data 
loggers  while  monitoring  the  performance  of 
solar  heated  buildings,    and  because  we  had 
an   extra  one  on  hand.     The  TRS-80  was 
connected  to  an   interface  which  incorporated 


an  8-channel,   8-bit  a/d  converter,    an  8-bit 
digital   output  port  and  an  8-bit   digital 
input  port.     All  ports  can  be  controlled  by 
commands   from  BASIC.     Four  bits  of  the  output 
port  were  used  to   control  the  rotation  of  the 
altitude   and  azimuth  motors  (two  bits   each). 
Both  the  azimuth  and  altitude  motors  were 
fitted  with  a  limit  switch  which  was  actuated 
when  the   sensor  was   in  the  "home"  position. 
Two  bits  of  the  input  port  were  used  to  sense 
these  switches.     A  real  time  clock  based  on 
a  quartz  oscillator  was  monitored  through 
another  bit  on  the  input  port.     This  input 
allowed  the  system  to  keep  track  of  time. 

Three  of  the  eight  a/d  converter  inputs  were 
used  for  this  experiment.     The   first   input 
was  used  to  digitize  the  signal  from  the 
scanning  Li-Cor  radiometer.     The  second  ana- 
log input  monitored  the  output  of  an  Eppley 
PSP  pyranometer  mounted  at   a  60°  tilt   facing 


Fig.  3. 
PSP. 


Scanning  radiometer  and  reference 
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due  south.     The  third  analog  ohaimal  moni- 
tored the  output  of  an  Analog  Devices  AI)590 
temperature  transducer.     Signal   conditioninR 
circuits  «ere  fatrioatsd  for  each  of  ihe 
transducers  which  amplified  and/or  offset 
the  transducer  signal  to  utilize  the  full 
5  V  range  of  the  a/D  converter.   ..The  com- 
pleted scanning  radiometer  and  the  reference 
Eppley  PSP  axe  shown  in  Figure  3.     The  com- 
puter,  interface,   eteppijig  motor  driver  and 
cassette  data  storage  are  shorn  in  Figure  4. 


Computer  and  interface 


''•0     ACCURACY  OF  SYRTFUl 

The  stepping  motors  used  in  this  system  are 
relatively  low  cost   and  are  manufactured  try- 
North  American  Phillips.     They   are  designed 
to  move   7.5°  or   15°  per  step.     The  nominal 
accuracy  of  these  motors   is   %  or  0.75°.     The 
outputs  of  the  solar  transducers  wore  ampli- 
iled  using  a  two-stage  operational   amplifier 
circuit  having  a  gain  of  about   370.     The 
temperature  sensitivity  of  the   ainplifier  cir- 
cuits and  the  analog  digital  converter  were 
itl  tt  ^  ?PPlyi"S  a  fixed  input   and  observ- 
ing the  output  while  the  temperature  of  the 
!  1?;?^  ""^^  assembly  was   forced  through 

Ih^Jt^rtt^'  difference.  These  tests  ^ 
showed  that  the  outputs  were  within  +O.5  bit 
for  a  temperature  change  of  30°  C.  fhe  ma^i- 
mum^reso  ution  of  the  8-bit  converter  i.  one 
P^  in  256,  or  about  0.%.  since  the  elec- 
tronics are  normally  located  within  a  "hJL" 


condition   space,   this  change   in  sensitivity 
was   judged  negligible. 

The   indexing  accuracy  of  the  altitude  motor 
was    .ested  on  a  clear  day  when  the   solar 
radiation  was  practically  constant   over  small 
time  intervals.     The  system  was  programmed 
to  scan  through  13  altitude  positions  and 
then  return  without  changing  azimuth.     The 

The^.^^r*'f  *r'  '"^^^^^'^  ^^   «^=h  position. 
The  results  of  this  test  showed  that  9%  of 

the  readings  of  solar  radiation  were  within 
■£.5  bit.     The  results  were  the  same  for 
^w!ri  ^fT"  °^  ^""Uth,   implying  that  the 
altitude  indexing  errors  were  within  the 
resolution  of  the  A^  converter. 

When  the  solaj-  transducer  is  at  an  incident 
angle  of  75°,   an  error  of  +10  „ould  produce 
a  solar  radiation  variation  of  +3  bits.     The 
measured  experimental  variations  of  +0.5  bit 
agree  closely  with  the  published  indexing 
accuracy  of  the  motor,  which  is  +O.40.     X 
similar  test  on  the  azimuth  motor  showed  a 
reproducibility  of  +1  bit.     This  value 
corresponds  to  the  nominal  accuracy  of  the 
motor,  which  is  0.75°.     The  combined  indexing 
errors  appear  less  than  +1°. 

The  overall  electrical  noise  in  the  assembled 
Bystani  was  tested  by  prograjnming  the  trans- 
ducer so  that  it  remained  at  a  fixed  angle 
on  a  clear  day  when  the  solar  radiation  was 
nominally  constant.     The  output  variation 
01  the  Bolar  transducer  as  indicated  by  the 
entire  system  was  constant  within  10.75  bits. 

The  test  results  show  that  errors  in  reading 
the  solar  radiation  with  this  apparatus  due 
to  errors  in  altitude,   azimuth,  noise  and 
temperature  drift  are,   on  the  average,  within 
7,^rtr  :   J^t  a  clear  day,   the  solar  radiation 
indicated  ly  the  radiometer  when  it  is  point- 
ed directly  at  the  sun  is  nearly  200  bits. 
These  accuracy  tests  indicate  that  errors 
due  to  all  sources  are  of  the  order  of  ^'jl,. 

The  references  cfuoted  previously   show  that 
the  silicon  cell   radiometer  is  typically 
within  a  few  percent  of  the  precision  pyra- 
nometer  and  often  within   1^  of  the  precision 
K^ranometer.     The  total  errors  incurred  in 
this  experimental  apparatus  are,  therefore, 
within  the  uncertainty  band  of  the  Bilioon 
ceil   radiometer. 


5.0     OFERATIOK   AKD  TEST  PROnKTVIRP 

The   scan   routine   is  depicted  in  Figure   5.     A 
single  scan  begins  with  the  axis  of  the 
transducer  horizontal  and  pointed  due   east. 
The  transducer  sequentially  increments   I50 
in  altitude  until   it   is  horizontal   and  facing 
due  west       At  this  point,   the  azimuth  incre- 
ments  150  towards  the  north  and  the  trans- 
ducer returns  in  15°  increments  until  the 
axis  is  horizontal  and  facing  south  of  due 
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Pig,  5.   Geometry  of  tjcan  routine, 

east  by  15°.  This  procedure  continues  until 
the  transducer  has  pointed  at  all  directions 
in  the  upper  hemisphere  on   15°  increments. 
The  transducer  then  returns  to  the  home  posi- 
tion of  due  east.  The  entire  scan  is  com- 
pleted in  less  than  30  seconds.  (The  scanner 
can  be  driven  at  over  2x  this  rate  but  the 
indexing  errors  will  increase  by  ■]%.) 

The  transducer  goes  through  a  horizontal 
position  on  each  180°  scan  of  the  azimuth 
angles.   If  the  solar  radiation  is  uniform 
throughout  the  period  of  the  scan,  then  the 
horizontal  readings  of  each  azimuth  scan 
will  agree  to  within  +1  bit.  The  multiple 
readings  of  the  horizontal  insolation  allow 
a  continuous  accuracy  check  on  the  indexing 
and  conversion  of  the  apparatus.  The  alti- 
tude and  azimuth  limit  switches  are  checked 
during  each  scan.  If  an  indexing  error  is 
detected,  the  data  from  that  scan  is  ignored 
and  the  sensor  is  automatically  repositioned 
oy  the  computer. 

At  the  end  of  each  scan,  the  temperature  and 
the  output  of  the  fixed,  Eppley  PSP  is  read. 
The  PSP  is  mounted  at  an  angle  of  60°  up 
from  horizontal,  facing  true  south.  This 
angle  corresponds  to  one  of  the  156  positions 
of  the  scanning  Li-Cor  radiometer.  Compari- 
son of  the  outputs  of  the  two  instruments  at 
this  angle  allow  a  further  check  on  the 
accuracy  and  validity  of  the  data.   ■ 

At  the  end  of  the  hour,  the  average  radiation 
at  each  of  the  156  scanning  positions  is  cal- 
culated. The  average  reading  of  the  fixed 
PSP  pyranometer  and  the  ambient  temperature 
are  also  calculated.  These  results  are  saved 
on  cassette  tape  along  with  the  Julian  day 
number  and  the  hour  of  the  day.  The  total 
data  base  each  hour  consists  of  I60  8-bit 
l)ytes.  The  data  tapes  generated  by  the 
Radio  Shack  TRS-80  microcomputer  are  subse- 
quently read  by  a  North  Star  Horizon  micro- 
computer and  the  data  transferred  to  floprar 
diskette.  These  data  files  are  stored  as 
0-bit  bytes  for  subsequent  processing. 

The  igystem  is  programmed  to  scan  and  save 
data  for  all  daylight  hours.  At  night,  the 
system  is  put  into  a  slower  scanning  mode, 
and  no  data  is  saved.  This  keeps  the  motors 
warm  and  prevents  snow  from  aocuimilating  on 
the  Li-Cor,  Figure  6. 


Fig.  6.  Effect  of  snow. 

6.0  SAMPLE  DATA 

A  video  display  of  the  raw  data  resulting 
from  one  scan  of  the  hemisphere  is  shown  in 
Figure  7.  The  matrix  of  hourly  average  data 


Fig,  7,  Display  of  output  during  a  single 
scan  on  a  cloudy  day. 

near  noon  on  February  24  is  shown  in  Table  1. 
This  day  was  clear  and  the  ground  was  covered 
with  snow.  The  format  of  this  Table  corre- 
sponds to  the  sequential  azimuth  and  altitude 
indexing  of  the  sensor.  This  output  is 
somewhat  "warped"  in  space  as  compared  to  the 
geometiy  of  the  readings  themselves.  The  sun 
was  almost  due  south  during  this  hour  so  the 
values  are  nearly  symmetric  about  the  0° 
azimuth  row. 

Table  2  shows  the  data  from  Table  1  normal- 
ized to  the  radiation  measured  at  a  60° 
south-facing  tilt.  During  the  past  five 
years,  we  have  measured  insolation  on  a  60° 
south-facing  tilt  at  30  locations  in  Montana 
(over  45,000  data  days),  (10).  He  hope  that 
data  from  this  scanning  experiment  will  help 
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TABLE  1 
lUTRlX  OF  HOURLY  ATERAOE  DATA  OH  A  CLKAH  DAY  HEAH  SOLAH  HOCHI,  kH-o"^ 

Altitude, deg. 

0    15    30    45    60    75    90   105   120   135   150   165   180 


-90 

0.40 

0.47 

0.54 

0.59 

0.63 

0,64 

0.64 

0.61 

0.56 

0.48 

0.40 

0.32 

0.31 

-75 

0.59 

0.66 

0.71 

0.73 

0.72 

0,69 

0.64 

0.56 

0.45 

0.34 

0.24 

0.26 

0.31 

-60 

0.76 

0.84 

0.88 

0.87 

0,82 

0.75 

0.64 

0.52 

0.36 

0.21 

0.20 

0.26 

0.31 

-45 

0.95 

1.01 

1.02 

0.98 

0.90 

0.78 

0.64 

0.47 

0.29 

0.15 

0.20 

0.26 

0.31 

i 

-30 

1.08 

1.13 

1.13 

1.07 

0.97 

0.84 

0.65 

0.44 

0.24 

0.15 

0.20 

0.26 

0.31 

■a 

-15 

1.16 

1.20 

1.19 

1,12 

1.00 

0.84 

0.64 

0.42 

0.19 

0.15 

0.20 

0,26 

0.31 

(S)   0 

1.18 

1,22 

1.21 

1.14 

1.02 

0.85 

0.65 

0.42 

0.19 

0.15 

0.19 

0.25 

0.30 

1 

15 

1.15 

1.19 

1.18 

1.11 

1.00 

0.84 

0.64 

o.ia 

0.19 

0.14 

0.19 

0.25 

0,30 

•H 

30 

1.07 

1.12 

1.11 

1.07 

0.96 

0.82 

0.65 

0.46 

0.26 

0.14 

0.19 

0.24 

0.30 

■31 

45 

0.95 

1.00 

1.01 

0.98 

0.90 

0.78 

0.65 

0.48 

0..30 

0.15 

0.19 

0,24 

0.30 

60 

0.78 

0.84 

0.87 

0.86 

0.82 

0.75 

0.65 

0.53 

0.39 

0.25 

0.19 

0.24 

0.29 

75 

0.61 

0.67 

0.72 

0.74 

0.73 

0.70 

0.65 

0.57 

0.47 

0.36 

0.26 

0,24 

0.29 

TABLE  2 
HOURLY  DATA  NORMALIZED  TO  60°  TILT  (30°  altitude)  AND  0°AZIMUTH,  M-m" 

Altitude, deg. 
0    15    30    45    60    75    90   105   120   1-?S   ISO   ifts 


-90 

0.33 

0.39 

0.45 

0.49 

0.52 

0.53 

0.53 

0.50 

0.46 

0.40 

0.33 

0.26 

0.26 

-75 

0.49 

0.55 

0.59 

0.60 

0.60 

0.57 

0.53 

0.46 

0.37 

0.28 

0.20 

0.21 

0.26 

-60 

0.64 

0.69 

0.73 

0.72 

0.68 

0.62 

0.53 

0.43 

0.30 

0.17 

0.17 

0.21 

0.26 

-45 

0.79 

0.83 

0.84 

0.31 

0.74 

0.64 

0.53 

0.39 

0.24 

0.12 

0.17 

0,21 

0.26 

a 

-30 

0.89 

0.93 

0.93 

0.88 

0.81 

0.69 

0.54 

0.36 

0.20 

0.12 

0.17 

0.21 

0.26 

Ti 

-15 

0.96 

0.99 

.?-?^i 

0.93 

0.83 

0.69 

0.53 

0.35 

0.16 

0.12 

0.17 

0.21 

0.26 

(S)  0 

0.98 

1.01 

i.OOl 

0.94 

0.84 

0.70 

0.54 

0.35 

0.16 

0.12 

0.16 

0.21 

0,25 

i 

15 

0,95 

0.98 

^nw 

0.92 

0.83 

0.69 

0.53 

0.35 

0.16 

0.12 

0.16 

0.21 

0.25 

•H 
t4 

30 

0.88 

0.93 

0.92 

0.38 

0.79 

0.68 

0.54 

0.38 

0.21 

0.12 

0.16 

0.20 

0.25 

«IJ 

45 

0.79 

0.83 

0.83 

0.81 

0.74 

0.64 

0.54 

0.40 

0.25 

0.12 

0.16 

0.20 

0.25 

60 

0.64 

0.70 

0.72 

0.71 

0.68 

0,62 

0.54 

0.44 

0.32 

0.21 

0.16 

0.20 

0.24 

75 

0.50 

0.55 

0.60 

0.61 

0.60 

0.58 

0.54 

0.47 

0.39 

0.30 

0.21 

0.20 

0.24 

establish  empirical  tilt  correlation  factors. 
These  factors  would  greatly  lucreaEe  the 
utility  of  the  existing  60°  tilted  data  base. 

As  an  aid  toward  visualizing  the  data,  models 
have  been  constructed  depicting  the  three- 
dimensional  surface  of  the  solar  radiation. 
The  photograph  in  Figure  8  shows  a  model  of 
the  solar  radiation  in  the  upper  hemisphere 


at  noon  on  a  clear  day  (data  from  Table  1). 
The  model  clearly  shows  a  large  "bulge"  in 
the  direction  of  the  sun.  The  intensity  of 
the  solar  radiation  drops  off  rapidly  for 
northerly  azimuths.  The  radiation  reaches 
a  minimum  when  the  transducer  is  pointed 
approximately  45°  above  the  horizon  and  due 
north.  As  the  axis  of  the  transducer 
approaches  horizontal,  additional  radiation 
is  recorded  due  to  reflections  from  the  snow. 
This  reflection  produces  a  "foot"  on  the 
model  when  the  transducer  points  north. 

Figure  9  shows  both  the  noon  model  and  a 
morning  model  (9i30  a.m.).  The  average 
-intensity  of  the  radiation  is  less  for  the 


Fig.  8.  Model  of  solar  radiation  in  upper 
hemisphere  at  noon. 


Pig.  9.  Morning  and  noon  models. 
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morning  model  but  the  shape  is  similar.  This 
model  also  shows  the  "foot"  due  to  hack- 
reflected  radiation  from  the  snow  covered 
surface. 


7.0  CONCLUSIONS 

(1)  The  results  presented  in  this  paper  are 
preliminary.  Nearly  two  months  of  satis- 
factory hourly  data  are  currently  on  file. 

(2)  At  the  writing  of  this  paper,  the  hugs 
in  the  apparatus  have  been  resolved  and  data 
is  being  logged  on  a  regular  basis.  Data 
will  he  collected  for  another  three  months 
under  the  present  contract. 

(3)  Future  work  will  include  statistical 
analysis  of  the  relationships  between  solar 
radiation  on  the  various  orientations  as  well 
as  comparison  with  existing  empirical  rela- 
tionships. 
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fhit  pi*er  describes  the  operation  and  reGulta 
of  a  current  state-supported  program,  Solar 
Insolation  Maasureraent,  Vnntaaa  (SIMM).     The 
goal  of  this  program  is  tc  meamrre  and  (tecru- 
■■ot  the  solar  insolation  available  at  30 
locations  withiE  the  stats  of  Montana.     To 
data,  oncer  25,000  solar  insolation  data-days 
iunra  bean  collected,     A  Biuanary  of  most  of 
this,  data  is  available  to  the  public  in  the  .    ■ 
iiontana  Solar  Data  Mantial ■  Fig.      1.     (l) 


montana 


solar    data 


SSsr=V^^^^ 


manual 


!mA 


Tig.      1,     Cover  of  Wontai^.t  Solar  Data  Manual. 
illiiatrs.tins  solar  insolaition  records  ( in— 
olndlng  eclipse  on  Feb.   16,   1979). 


Early  in  t977  the  SIMM  program  (Solar  Insol- 
ation Meaaurement,  Montana)  was  initiated  by 
PowUcas  Engineering  under  a  grant  frrmi  the 
Altamativo  Reaawable  Energy  Sources  Program 
of  the  Montana  Department  of  Natural  Re- 
Bources  and  Conservation.     Twenty  instruments 
to  measure  solar  radiation  ware  sited  across 
Montana  during  1977  and  an  additional  ten  in- 
Btrumenta  ware  sited  early,  in  1978.     In  add- 
ition to  the  thirty  stations  having  contin- 
■atma  insolation  recorders  there  are  currently 
twenty-five  stations  which  report  the  readingo 
of  manual  inatrtuaantB,     Fig.  2. 


a     Coatlauouo  Data  Recording  Station 

O     Manual  Data  Station 

Pig.     2,     Montana  Solar  Insolation  Keasura- 
mant  ^etwoiic 

All  of  the  instruaents  are  located  in  public 
hl^  schools.     At  each  school  there  is  a 
volunteer,  usually  a  science  teacher,  who  is 
reoponmible  for  the  instrument.     The  volunteer 
periodically  replaces  the  chart  paper  in  the 
racordar  and  nails  the  completed  strip  charts 
to  our  office  for  prooeoslne.     The  total  time 
required  for  the  volunteer  to  do  these  tasks 
usually  suaotmts  to  lens  than  fiv®  hours  each 
year. 

The  solar  insolation  strip  chart  recorders 
are  often  located  in  a  laboratory,  office,  or 
olassrooa  where  they  can  be  seen  by  the  stu- 
dents.    Sone  teachers  make  specific  use  of  the 
inaolation  records  in  their  classes. 
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ill  hl^  sebs«le  portlslpa-tlB^  la  the  SDtH 
progr«B  reoeiT®  a  copy  of  Tha  Handbook  of 
Solar  Energy  Erperlffien'tB  and  &  mauu&l   oolar 
iBfflolation  metar.     Tins  Handbook  was  vrittt  ■ 
liy  tba  author  of  this  paper  as  an  integral 
part  of  tha  Insolation  msasureaant  prograa. 
Tha  80  pa^  Handbook  contiuinB  a  rovisM  of 
Mlar  taaica  and  tan  outlined  ezperiaents 
>«nglTig  froH  Insolation  oeoBurement  to  test- 
ins  *  solar  collector. 


8.    tssoLiTias  an  a  tilted  5uh?ace 

■•taorologlcaJ.  iaooliitlon  data  ia  naimlly 
recorded  on  a  horliotital  surface  due,   in  p&rt, 
te  the  fact  that  the  solar  colloctor  of  in- 
tarart  to  metaorologistis  ia  the  horizontal 
,  nrfae*  of  the  earth.     For  northcm  U.S.   Itt- 
itndta  the  solar  designer  wall  normally  tilt 
tts  collector  plana  between  45°  and  90°  up 
Srvm  the  horisontal  to  increase  the  aryerage* 
radiant  flus  donaity  intercepted  by  the  coll- 
ector.   Procedures  exist  for  estimating  the 
insolation  on  a  tilted  surface  given  the 
.BBaoured  insolation  on  a  horiiontal  surface. 
Sowever,  the  net  uncertaintiea  of  these 
,_  estimates  io  northern  U.S.   latitudes  during 
■  the  urinter  will  be  around  *•  20^.     Some  un- 
cartaintiSB  arise  from  the  pyTanoooter  errors 
«t  low  solar  altitudes.     Additional  'jooer- 
taintles  eilBt  duo  to  reflected  radiation 
fron  snow-covered  ground  which  will  ba  In- 
tsrctpted  by  steep  solar  collectors  but  will 
aet  be  intercepted  by  a  horizontal  radiation 
transducer. 

81^;^  the  goal  of  the  SD3?  prograo  ia  to 
dwrolop  c  polar  insolation  data  base  for  de- 
■ijping  solar  conversion  devices  we  tneaaure 
insolation  on  a  gurface  tilted  60°  np  from 
horizontal  and   facing  due   routh.     The  trans- 
daeer  hence  becomes  a  model  01   a  typical  so- 
lar collector  which  might  be  uBed  at  this 
latitude  to  collect  heat  ia  the  winter.     Th« 
resulting  60°  Bolar  radiation  data  can  often 
lia  used  directijf  or  cjia  ba  tranefonoed  to 
aaigKboring  tilt  angles  with  relatively  naall 
inexaaseB  in  uncertainty.     Tha  effects  of 
Boow  sad  dirt  on  a  collector  are  approiii- 
aatslj-  modeled  by  this  transducer  orientation 
■nd  srs  an  inta^Ta  pairt  of  the  iasolation 
data. 


3.     CTSTRPKEOTATIOir 

Shs  solar  insolation  is  isensed  by  a  silicon 
c«ll  Bounted  in  an  Ban  square  aluminum  plate 
■  and  covered  with  glasB,     The  transducer  is 
attached  to  a.  bracket  outuide  the  building 
Bad  »  wire  runs  to  a  strip  chart  recorder 
located  inside  the  building.     Ve  fabricated 
oil  th»   transducers  used  in  this  program. 

Sllioon  cell  charactoristlcB  relating  to 
hoLbt  flni  oeaBuranant  jire  wall  docura«ntod  in 
t&a  literature  and  aovoral  coosaarcial  silicon 


P«ll  iastruaantB  are  wrailabls  (2).     Charact- 
•ristles  OBUslag  departuraa  from  pyranometar 
roadinga  are  dus  to   spectral  rsBponse,   temp- 
exaturo  ciajgai  and  coBine  rsBponae.     The 
Bp«otr8l  re.Tponsa  will  cauBO  ths  silicon  cell 
to  read  ooajrurably  higher  on  cloudy  days. 
Ths  instmcenis  uaed  in  this  program,   for  es.-  ■ 
•Bpls,  will  raad  about  0.01 5  Wh/o^  high  .under 
'  conditions  ot~  'iCff,  beas  attenuation  due  to 
oloudffl.     The  effect  of  this  spectral   error  on 
the  ■onthly  average,   daily  total  energy  dens- 
ity is  typically  reducaa  to  less  than  3^. 
TM«  r«dac»d  error  is  daa  to  the  relatively 
loM  eosrgy  density  on  e^loud  days.     Cosine 
•rrom  are  minimized  dae  to  the  tilted  orient- 
ation. 

Thffl  olisotrieal  circuit  tuied  is  sensitive  to 
changes  in  lead  resistaacu.     This  factor  has 
probably  been  a  primary  raource  of  error  in 
this  maaauretaent  prograsa,     Raaistanoe  vari-. 
ationa  are  usually  dus  to  improperly  soldered 
lead  wire  connections.     The  gradually  in- 
creasing contact  resiBtsBice  leads  to  a  grad- 
ual reduction   in  sensititrity  of  the  instru- 
ment.     In  a  few  caaes  this  problem  has  been 
responsible  for  temporary  errors  of  10  and 
.  20^.     He  are-  currsmtly  iaspectijig  andsarvio-  ■  . 
izig  ths  instruments  at  mil  30  stations. 


4.     I&TA  REPCCTION 

Strip  chart  data  is  processed  by  tracing  ths 
curves  on  a  digitizer  paaten  connected  to  a 
microconpirter,  Pig.   3.     The  resolution  of  the 
platen  is  0,1  mo  and  each  curve  is  represented 
'bj  200  to  400  coordinatas.     The  reading  error 
du®  to  resolution  is  less  than.1^.     The,  prim- 
ary oourco  of  data  redaction  error  is  the 
repeatability  of  the  curve  tracing  itself. 
On  clear  days,  where  th»  curve  is  a  well  de- 
fined line,   the  reproducibility  averages  1%. 
On  cloudy  dayti  the  reproducibility  can  reach 
5  to  20%.     This  is  due  to  the  fast  response 
■of  tha  silicon  cell  which  produces  a  very 
rough,   jagged  curve  on  ooms  cloudy  days. 
Fig.   U 


Pig.   3  Computer  |ayctea  xmrni  to  read  strip-charts 
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Itarlnc  ■til*  1b«*  r«8*  r«produci1jili-ty  taste 
h«ra  baan  run  on  over  800  data  da^s  (^  of 
OOP  total  data  base) .     Th«  roproduoibility  of 
tha  Bonthly  average,   daily  total  energf  dene- 
itf  »Tr«rag«»  better  than  -  2%.     Theoa  montlily 
Vkluaa  often  reproduce  within  1?5  duo  to  th® 
eoopenaation  of  ];>ositi7a  and  negativQ  errors 
whan  »varag«d« 

Tha  procaaaod  data  from  each  Aar  consists  of 
dv  nuahar,   total  energy  daneity,   day  lenfrth, 
aooa  maxiraum  intensity  ( if  day  wao  clear) ; 
in'tenaity  Bpactrura,   half-iour  average  aolar 
Intmaitles,   and  a  faiilt  coda  (specifying 
loat  or  suspect  data) .     The  total  energy  is 
atored  in  2  (3  hit)   bytes  and  thus  has  a 
Baxlnmn  resoiution  of  1   part  in  512.     All       ^ 
other  non-integar  cjuantities  are  rounded  to 
intogora  and  stored  in  a  single  byto  in  the 
fila  (oaximin  resolution  1  part  in  256).     The 
m^ninvi  X0S8  of  precision  in  this  data  due  to 
rounding  off  is  0. 5^. 

Tha  data  are  stored  on  hard-«®ctorad,   5i  inch 
■a^atlo  diskettes  using  tha  Morth  Star  Disc 
'  Oparating  System     (a  popular  miei^ocoinpater.  ■  '•  ■ 
fomat).     Each  single-<iensity  diakatts  con- 
tains 3  statlon-^ears  of  solar  data.     This 
dftta  has  baen  used  by  PcwUces  Engineering  aa 
input  to  solar  heating  simulation  models. 
Coplaa  of  thesa  diskettes  can  ba  supplied  at 
B«)fif;t»ifil   coat  for  users  wanting  to  impleaant 
thair  own  aimolation  prograas. 


5,     DATA  ACetmACT.   CLEAR  BAY  METHOD 

All  traaaducers  are  calibrated  bofore  install-        Pig- 
ation  to  establish  a  nominal  calibration 
factor.     Tha  reference  instruments  haire  in- 
cluded Kipp  pyrsnometers,   Li-Cor  radio- 
aetsrs  and  one  of  our  transducers  which  was 
calibrated  by  KOAA.     This  provides  an  initial 
rafarence  point  but  does  not  allow  for  tha 
effects  of  drift  and  aging,  which'  may  occur 
while  the  instrunent  is  on  site.     For  long 
tana  control  of  accuracy  we  make  oaa  of  a 
oXaar  day  technique  similar  to  the  methoda 
used  to  adjust  the  SOIXE?  data  base  ( 3) . 

Tha  aaaanee  s'f  our  clear  day  prooadurea  is  to 
eoopara  tha  maajsurad  radiant  fluz  at  noon  to 
a  value  predicted  by  a  theoretical  clear  day 
■odel.    He  have  studied  two  clear  day  tnodela. 
The  first  model  was  provided  by  Douglas  Hoyt 
(4)    and  the  second  model  was  baaed  on  the 
ASHRAE  otpiationa  with  a  claamess  modification 
to  accannt  for  altitude. 

^  part  of  our  data  reduction  procedures  we 
pick  off  tho  noon  mg|-riTmT?n  OQ  all  clear  days 
and  this  value  is  stored  in  the  data  file. 
Before  suimnarizing  the  annual  data  wo  have  a 
program  that  scans  all  clear  days  and  computaa 
a  theoretical  value  for  that   «amo  day.     If 
there  is  a  discrepancy  of  mora  than  a  few  per- 
cent,  the  eiperimental  data  is  scaled  to  agree 
with  th«i  theoretical  value.     Figure  4  shows  a 


oMQarleon  of  thsorotloal  and  azperlaantal 
clear  day  values  for  two  stations  in  1978. 
Piguro  4  (b)   shows  data  that  was  subsequently 
scaled  for  a  batter  fit.     Ha  feel  that  this 
prooedura  Inanros  that  our  data  base  doas  not 
contain  significant  or  large  errors.     Aitar 
scaling  most  of  tha  stations  match  tha  theor- 
etical clear  dsy  values  to  withing  t.  2^  or 
better. 
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(a)  Qt.  Palls  data  for  1978 
shoMing  two  clear-day  models  with  data 


Bodified  ASHRAE  modal 
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Cb)  Glendive  data  for  1978 
illumtratlng  use  of  seals  factor 


4    Compoxison  of  clear  day  experimental 
azd  theoretical  values 


The  solar  insolation  data  from  this  program 
has  been  suaiBarized  in  the  Montana  Solar  Data 
Manual.     The  first   edition  of  this  manual  was 
distributed  to  libraries,   extension  services, 
anginesra,   and  architects  in  the  aummer  of 
1979.     An  updated  edition  will  bo  distributed 
in  1980.     The  now  manual  includes  thrae  years 
of  solar  insolation  data  in  tables  of  weakly 
and  monthly  averages  of  daily  totals.     Also 
included  are  tables  of  the   frequency  of  occur- 
anoo  of  oloudy  days  and  a  reproduction  of 
oliiaatological  data  from  NHS   records.     Tha 
mgr;<ni1    includes  a  discussion  of  solar  heating 
design  mtthodffl,    economic  analysis,    a  list  of 
raferencea,    conversion  tables  and  a  discussion 
of  the  accuracy-  of  tha  solar  data.     The  manual 
containB  about  200  pages  of  information  and  is 
availahlo  from  the  state  agency  at  no  cost. 

Tha   average   intensity  speetroa  is  also  pre- 
sented aa  a  percentage  of  the  daily  total 
energy  density  which  occurod  at  radiation   in- 
tensities above  12  threshold  values  ranging 
from  0.1    to   1.2  MJ/m^, 


nqo 


datft  frwa  th«  Montana  Solar  Data  Manual   for 
H»l«l*f  Kontana  d-jxing  1978-     The  first  coi- 
nan  llata  tbo  Btarting  date  for  the  weak  and 
th«  Bscond  colunm   listo   the   avera^   daily 
total  «n«rsr  denaity,   kW-hr/m'--<iay9   duxiYig 
that  week.     Tha  next  twelve   coluxao   list  the 
vrara^  apectrua  during  the  weeic,     Sioiilarly, 
»or«  compact  tables   listing  monthly  averages 
9T9  also  included  in  the  Manual,      Parentheses 
around  the  daily  avera^s  denote  that  at  least 
oa«  d«y  of  data  vas  ndssing  d:iring  the  avor- 
ft^4j  period* 
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?ig.   5    Sasupla  solar  dsta  fS-oa  ■aaual 

Til*  cnrrs3a"t  Manual  coniaojiB   awlar  data  frxjiu 
the  fiTBt  ^hreffl  yaam  of  •the  aaaEurtamsat  paro— 
gris,    1977}    1978  acd  1979.     TJiess  moaBuromBa'tB 
ar«  'baing  conticuad  througii  the  aaxt  year  and 
thB  Ds-ta  Kau'J&l  will  ba  UF<isted  at  that  tiaa. 
If  fonda  ar®  a-ffailatle,    this  prograa  will  bo 
eontiausd  to  a.  point  whsre  long-tana  areras* 
iasoLation  data  can  be  raliably  docimanted. 
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